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ABSTRACT 


A  Time-Dependent  Model  for 
the  Low-Latitude  Ionosphere 

by 

Gary  D.  Wells,  Master  of  Science 
Utah  State  University,  1988 

Major  Professor:  Dr.  Robert  W.  Schunk 
Department:  Soil  Science  and  Biometeorology 

V, 

The  theoretical  development  of  a  time-dependent,  self-consistent  model  of 
the  low-latitude  ionosphere  between  200  and  3000  km  is  presented.  This  model 
solves  the  coupled  continuity-momentum  and  energy  equations  for  an  ionosphere 
composed  of  two  ion  species  (O"^  and  H"*")  and  electrons  in  a  quasi-simultaneous 
manner.  Account  is  taken  of  effects  due  to  the  curvature  of  the  field  lines  and  the 
inclination  of  the  field  line  to  the  vertical.  This  model  allows  for  electrodynamic 
drift  of  the  plasma  during  the  day,  effects  due  to  the  expansion  and  contraction 
of  the  neutral  atmosphere,  and  ion-neutral  drag  effects,  as  well  as  the  transfer 
of  momentum  through  collisions.  The  energy  equation  allows  for  heating  or 
cooling  through  collisions,  conduction  parallel  to  B,  frictional  heating/cooling 
due  to  the  relative  motion  between  the  ions  and  the  neutrals,  and  effects  due  to 
adiabatic  expansion  and  contraction.  '  / 

We  have  significantly  improve'tf' modelling  capabilities  by  including  more  re- 

/■ 

alistic  neutral  winds  and  electric  fields.  We  have  also  allowed  for  the  offset  in 
the  geographic  and  geomagnetic  poles  by  mapping  the  neutral  winds  and  time 

r 

(both  organized  in  the  geographic  coordinate  system)  into  the  magnetic  refer¬ 
ence  frame.  By  mapping  the  neutral  winds  into  the  magnetic  coordinate  system. 


longitudinal  variations  in  the  density  and  temperature  of  the  F-region  are  intro¬ 
duced  in  a  consistent  manner.  The  net  result  is  a  model  that  has  applicability 
for  a  variety  of  local  times,  seasons,  longitudes,  and  altitudes. 

Results  from  this  model  are  shown  to  reproduce  the  important  physical  pro¬ 
cesses  that  are  unique  to  the  low-latitude  ionosphere,  such  as  the  depletion  of 
plasma  at  the  magnetic  equator  due  to  electrodynamic  drifts  and  diffusion,  and 
the  seasonal  anomaly. 

(110  pages) 
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CHAPTER  I 


INTRODUCTION 


The  Terrestrial  Equatorial  Ionosphere 

The  terrestrial  ionosphere  is  that  portion  of  the  earth’s  upper  atmosphere  in 
which  ions  and  free  electrons  exist  in  significant  numbers.  This  is  not  to  say  that 
ion  densities  exceed  the  neutral  densities  throughout  the  ionosphere  for  in  the 
lowest  regions  of  the  ionosphere,  the  neutral  concentration  is  still  many  orders 
of  magnitude  greater  than  the  ion  concentration.  In  fact,  ion  densities  generally 
do  not  exceed  neutral  densities  below  1000  km  for  moderate  solar  activity. 

The  upper  extent  of  the  ionosphere  is  often  referred  to  as  the  plcismasphere. 
The  plasmasphere  is  a  torus-shaped  region  which  encircles  the  earth  at  low  and 
middle  latitudes.  It  is  composed  of  a  cold,  relatively  dense  plasma  of  ionospheric 
origin  which  corotates  with  the  earth.  In  addition  to  corotating,  the  plasma  can 
flow  between  conjugate  hemispheres  as  depicted  in  Figure  1.  The  plasmasphere 
has  a  maximum  vertical  extent  of  4  earth  radii  (R^)  at  the  equator.  At  even 
greater  heights,  a  transition  from  plasmaspheric  to  the  magnetospheric  plasma 
occurs.  This  transition  is  called  the  plasmapause  and  is  marked  by  a  large  de¬ 
crease  in  electron  density  [Schunk,  1983).  Young  et  ai,  [1979,  1980]  and  Richards 
et  al.  [1983]  have  developed  time-dependent  models  valid  for  plasmaspheric  fleld 
lines,  along  which  the  important  variations  are  in  the  vertical  direction.  This 
makes  this  region  somewhat  easier  to  model  than  the  lower  ionosphere,  where 
variations  in  the  horizontal  and  the  vertical  are  important  [Schunk,  1988). 

Our  region  of  interest  will  be  the  low-latitude  ionosphere  where  some  unique 
transport  effects  occur  due  to  the  near  horizontal  nature  of  the  field  lines.  The 
first  effect  is  due  to  thermospheric  winds  blowing  from  the  surmner  to  winter 
hemisphere.  These  winds  induce  a  drift  in  the  ionospheric  plasma;  and  since 
the  ions  and  electrons  are  bound  to  the  magnetic  field  lines,  this  drift  is  field- 
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Figure  1.  Schematic  representation  of  interhemispheric  plasma  flow  along  geo¬ 
magnetic  field  lines  crossing  the  equator  at  plasmaspheric  altitudes  for  solstice 
conditions  (from  Evans  and  Holt,  (l978j).  Reprinted  by  permission  of  D.  Reidel 
Publishing  Company. 

aligned  (This  effect  is  shown  schematically  in  Figure  2).  As  the  plasma  rises  in 
the  summer  hemisphere,  it  expands  and  cools  adiabatically;  and,  conversely,  as 
it  descends  in  the  winter  hemisphere,  it  is  compressed  and  warmed. 

Another  feature  is  the  equatorial  fountain  or  Appleton  anomaly  which  forms 
during  the  day  on  both  sides  of  the  magnetic  equator.  This  feature  consists  of 
two  peaks  in  the  latitudinal  profile  of  electron  density  centered  on  either  side  of 
the  magnetic  equator,  and  a  “trough”  in  the  electron  density  at  the  magnetic 
equator.  In  the  da>  time  .^-region,  thermospheric  winds  induce  ionosphei  ic  cur¬ 
rents,  which  induce  an  eastward  electric  field  in  the  F  region  (See  Figure  3  for 
the  approximate  altitudes  of  the  various  ionospheric  regions) .  The  combination 
of  this  eastwaurd  electric  field  and  a  northward  magnetic  field  produces  an  up¬ 
ward  drift  of  the  plasma  (ions  and  electrons)  at  the  magnetic  equator.  As  the 
plasma  is  being  lifted,  it  is  acted  on  by  gravity,  which  causes  it  to  drift  away 
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Figure  2.  A  schematic  representation  showing  the  vertical  drift  of  the  ionospheric 
plasma  (VF)  due  to  a  horizontal  neutral  wind  along  a  magnetic  meridian  (C/^). 
The  inclination  of  the  magnetic  field  to  the  horizontal  is  defined  as  I  (from 
Bittencourt  and  Sahai,  [1978]).  Reprinted  by  permission  of  D.  Reidel  Publishing 
Company. 
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Figure  3.  A  schematic  illustration  of  the  heights  and  densities  of  the  various 
ionospheric  regions  (from  Banks  et  al.,  [1976]).  Reprinted  by  permission  of  D. 
Reidel  Publishing  Company. 
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from  the  magnetic  equator  and  diffuse  down  the  magnetic  field  lines  forming  the 

0 

regions  of  enhanced  plasma  density  at  ±15  magnetic  latitude  and  depleting  the 
plasma  density  at  the  magnetic  equator.  This  process  is  shown  in  Figure  4. 
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Figure  4.  An  illustration  of  the  processes  leading  to  the  formation  of  the  Apple- 
ton  peaks  in  the  daytime  F-region  at  equatorial  latitudes  (from  Schunk,  [1983]). 
Reprinted  by  permission  of  D.  Reidel  Publishing  Compamy. 


Historical  Development 

Determination  of  the  important  factors  controlling  the  formation,  mainte¬ 
nance,  and  variation  of  the  F-region  and  the  topside  ionosphere  have  resulted 
in  the  development  of  a  number  of  theoretical  models  over  the  last  24  years. 
The  decade  of  the  1960’s  was  a  period  in  which  the  primary  focus  was  on  estab¬ 
lishing  those  processes  which  control  the  behavior  of  the  F-region.  Diffusion, 
while  an  important  process  in  this  region,  was  shown  to  be  unable  to  account 
for  the  formation  of  the  equatorial  anomaly,  the  maintenance  of  the  ionosphere 
at  night,  and  the  seasonal  anomaly  \Rishbeth  et  ai,  1963;  Rishbeth,  1964,  1966] 


even  when  the  diurnal  expansion  and  contraction  of  the  neutral  atmosphere  and 
the  ionosphere  were  taken  into  account  [Thomas,  1966;  Thomas  and  Venables, 
1967;  Rishbeth,  1964).  Their  results  did  show  that  the  peak  in  the  i^-region 
formed  at  the  altitude  where  diffusion  is  balanced  by  photochemical  loss,  that 
diffusion  was  the  dominant  process  above  the  peak,  and  that  photochem¬ 
istry  was  dominant  below  the  F2  peak  [Rishbeth  et  al.,  1963].  Variations  in  the 
temperature  were  found  to  affect  ion  densities  and  the  maximum  height  of  the 
F2-region  [Thomas  and  Venables,  1967];  however,  thermal  expansion  and  con¬ 
traction  of  the  neutral  atmosphere  has  little  effect  on  ion  densities  during  the 
day  or  night  [Schunk  and  Walker,  1973]. 

In  the  early  1960’s,  the  presence  of  an  electric  field  within  the  F-region  had  yet 
to  be  confirmed  -  but  was  thought  to  exist  -  [Bramley  and  Peart,  1965]  and  some 
of  the  diffusion  models  being  developed  assumed  the  magnetic  field  lines  were 
approximately  vertical.  Easier  (1964,  1967],  Sterling  et  al.  [1969],  and  Kendall 
[1962]  were  largely  responsible  for  developing  the  appropriate  equations  and 
diffusion  coefficients  allowing  for  calculations  along  magnetic  field  lines  where 
the  coordinate  system  (or  field  line)  was  moving  with  the  vertical  £  x  B  drift  of 
the  plasma.  To  an  observer  on  the  the  field  line,  this  would  mean  that  all  plasma 
motions  appeared  to  be  parallel  to  the  magnetic  field.  By  eissuming  a  geocentric 
dipole  in  which  the  earth’s  rotational  and  magnetic  axes  are  aligned  and  diffusion 
along  the  field  line,  it  was  demonstrated  that  even  relatively  small,  upward 
vertical  drifts  (10  -  20  meters/second)  at  the  equator  combined  with  diffusion 
down  the  field  lines  away  from  the  equator  could  account  for  the  formation  of 
the  equatorial  anomaly  during  the  day  [Easier,  1964;  Bramley  and  Peart,  1965]. 
The  diurnal  behavior  of  the  equatorial  anomaly  was  reproduced  by  assuming  a 
sinusoidal  variation  of  this  electric  field  [Baxter  and  Kendall,  1968]. 

The  processes  controlling  the  asymmetrical  distribution  of  the  equatorial  anom¬ 
aly  and  the  nighttime  maintenance  of  the  F-region  had  yet  to  be  defined.  Spec¬ 
ulation  concerning  the  nighttime,  low-latitude  F-region  involved  two  processes. 


The  first  considered  an  external  source  of  ionization  such  as  an  influx  of  en¬ 
ergetic  particles  or  an  influx  of  ions  from  the  protonosphere.  The  second 
proposed  a  change  in  the  rate  at  which  the  ions  recombine  either  by  a  change 
in  the  composition  of  the  nighttime  ionosphere  or  by  vertical  movement  of  the 
F-region  into  a  region  of  lower  atmospheric  density.  Computation  of  the  day¬ 
time  flux  of  H"*"  ions  to  the  protonosphere  showed  that  the  protonosphere  could 
not  provide  the  necessary  downward  flux  of  ions  to  maintain  the  F-region 
at  night.  In  addition,  rocket  measurements  indicated  an  upward  motion  of  the 
F-region  plasma  [Hanson  and  Patterson,  1964].  This  upward  motion  is  the  re¬ 
sult  of  an  equatorward,  meridional  neutral  wind  which  becaiise  of  drag  between 
the  neutrals  and  ions  forces  the  plasma  up  the  inclined  magnetic  field  lines  at 
low  latitudes  [Rishbeth  and  Garriott,  1969].  The  meridional  asymmetry  of  the 
equatorial  anomaly  was  effectively  modelled  by  the  inclusion  of  a  horizontal  neu¬ 
tral  wind  blowing  from  the  summer  to  winter  hemispheres  for  solstice  conditions 
into  models  containing  diffusion  and  electrodynamic  drifts  [Sterling  et  al.,  1969; 
Abur-Robb  and  Windle,  1969;  Bramley  and  Young,  1968],  The  longitudinal  vari¬ 
ation  of  the  equatorial  anomaly  can  be  obtained  by  allowing  for  the  offset  of 
the  magnetic  and  rotational  cixes.  This  offset  shifts  the  neutral  wind  pattern 
-  organized  in  geographic  coordinates  -  in  relation  to  the  magnetic  field  lines 
and  in  conjunction  with  a  longitudinal  variation  in  the  F-region  electric  field 
produces  longitudinal  variations  in  the  equatorial  anomaly  [Anderson,  1973a, 
19736;  Chan  and  Walker,  1984a,  19846). 

In  1975,  Schunk  [1975],  noting  that  the  basic  trzinsport  equations  used  in 
aeronomical  research  had  not  generally  been  applied  with  a  clear  recognition 
of  their  limitations,  developed  a  general  system  of  transport  equations  valid  for 
a  multi-component  gas  mixture.  This  system  of  equations  includes  the  conti¬ 
nuity,  momentum,  internal  energy,  pressure  tensor,  and  heat  flow  equations  for 
each  species  of  the  geis.  The  basic  problem  for  aeronomical  studies  is  choosing 
the  appropriate  set  of  equations  for  the  velocity  distribution  inherent  to  each 


7 

species  in  the  gas.  For  undisturbed  conditions,  a  Maxwellian  distribution  is 
appropriate,  but  the  presence  of  various  non-equilibrium  processes  often  causes 
deviations  from  the  MzocweUian  distribution.  Therefore,  before  the  moments 
of  the  Boltzmann’s  equations  can  be  obtained,  the  velocity  distribution  must 
be  specified.  Generally,  the  gas  mixture  has  been  assumed  to  be  Maxwellian 
and  the  fluid  equations  have  been  applied.  This  general  system  of  transport 
equations  is  developed  using  a  13-moment  approximation  and  is  valid  for  both 
the  neutral  atmosphere  and  the  ionosphere.  They  can  be  applied  to  collision 
dominated  and  collisionless  plasmas,  and  in  the  case  of  very  frequent  collisions, 
these  equations  give  the  usual  Navier-Stokes  expressions  for  the  viscous  stress 
tensor  and  the  heat  flow  vector. 


During  the  late  70’s  considerable  attention  was  given  to  the  problem  of  inter- 
hemispheric  flows  of  plasma  along  dipole  field  lines.  Time  dependent  and  steady 
state  solutions  were  obtained  from  the  continuity  and  momentum  equations  be¬ 
tween  /^-regions  in  the  conjugate  hemispheres.  Bailey  et  al.  [1977o]  studied 
the  interhemispheric  flows  of  0+  and  ions  for  sunspot  minimum  conditions 
at  equinox  at  L  =  3.  They  found  that  counterstreaming  between  the  and 
H"^  could  indeed  occur  during  the  night  at  solar  minimum,  and  to  account  for 
the  differences  between  their  calculated  values  and  observed  values,  they  con¬ 
cluded  electrodynamic  drift  was  important  at  the  time  the  observations  were 
made.  They  also  studied  the  nighttime,  field-aligned  fluxes  at  solar  maximum 
for  L  =  3  [Bailey  et  al.,  19776).  They  showed  that  the  O'''  flux  was  not  related 
to  the  flux  of  H'*'  to  or  from  the  protonosphere;  and,  hence,  nighttime  F-region 
densities  are  not  strongly  affected  by  fluxes  to  or  from  the  protonosphere. 


In  1979,  Bailey  and  Moffett  [1979]  extended  the  previous  calculations  to  in¬ 
clude  the  flux  of  He”^  ions  (in  addition  to  H"''  and  O''')  and  thermal  effects  on 
field-aligned  fluxes  by  performing  steady  state  calculations  of  densities,  tem¬ 
peratures,  and  flow  velocities  for  each  species.  The  ions  were  heated  through 
collisions  with  the  electrons,  the  neutrals,  and  each  other;  and  the  electrons  were 


heated  through  collisions  with  photoelectrons,  Jind  cooled  due  to  the  fine  struc¬ 
ture  excitation  of  neutral  O,  and  the  vibrational  and  rotational  excitation  of  O2 
and  Nj.  They  assumed  a  geocentric  dipole  field  line  with  an  equatorial  crossing 
altitude  of  930  km  and  a  north-to-south  neutral  wind.  They  concluded  that 
the  transequatorial  breeze  is  not  significantly  affected  by  elevated  plasma 
temperatures,  but  the  field-aligned  velocities  of  and  He"*"  were  strongly  de¬ 
pendent  on  temperature  gradients  and  the  O'*"  field-aligned  velocity.  Therefore, 
for  the  light  minor  ions  (H'*'  and  He"^)  ion  -  ion  drag  cannot  be  ignored.  It 
has  also  been  reported  that  the  £[■*■  flux  is  upward  throughout  the  day  and  not 
significantly  affected  by  the  H"^  tube  content;  whereas,  at  night  the  flux  can 
be  somewhat  variable.  The  daytime,  peaJc  density  in  the  F2  region  is  reduced 
by  about  15%  by  this  upward  flux  of  and  the  nighttime  peak  density  in  the 
F-region  is  maintained  by  the  neutral  winds  and  the  downward  flux  of  [Mur¬ 
phy  ef  al.,  1976].  In  a  study  of  counterstreaming  between  O"*"  and  H"*",  Young  et 
al.  [1979]  considered  the  case  of  a  dipole  field  tilted  with  respect  to  the  earth’s 
rotational  axis  under  steady  state  conditions  and  found  that  counterstreaming 
must  occur  under  these  conditions  as  long  as  the  neutral  and  ion  temperatures 
are  not  the  same. 

A  final  area  of  investigation  involves  the  effects  of  interhemispheric  transport, 
electrodynamic  drifts,  and  heating  by  photoelectrons  on  the  plasma  tempera¬ 
tures.  Early  studies  by  Bailey  et  al.  [1973,  1975]  involved  the  steady  state 
solution  of  the  energy  equation  with  densities  being  supplied  by  the  model  of 
Moffet  and  Hanson  [1973].  The  ions  were  assumed  to  be  heated  through  colli¬ 
sions  with  the  electrons,  cooled  by  collisions  with  the  neutrals,  and  to  have  the 
same  temperature  (Tc)+  =  T/^+).  The  electrons  were  heated  through  the  exci¬ 
tation  of  neutral  O,  O2,  N2,  and  the  interaction  with  photoelectrons  [Bailey  et 
al.,  1975].  These  solutions  were  valid  for  daytime  and  nighttime  during  equinox 
for  low  and  moderate  solar  activity.  Later  studies  [Bailey  and  Heelis,  1980;  Bai¬ 
ley,  1980;  Bailey,  1983;  Moffett  et  al.,  1986]  were  performed  with  self-consistent, 


time-dependent  continuity,  momentum,  and  energy  equations  with  an  empirical 
expression  for  the  heat  flow.  The  results  for  equinox,  solar  maximum  with  no 
interhemispheric  transport  and  no  electrodynamic  drift  above  Arecibo  (L  =  1.4) 
indicated  that  for  altitudes  between  120  km  and  2548  km  was  the  dominant 
ion  during  the  day,  but  was  comparable  to  at  night.  During  the  night  the  ion 
temperatures  are  approximately  the  same,  but  during  the  day  the  tempera¬ 
ture  can  exceed  the  O'*"  temperature  by  300°K  [Bailey,  1980].  For  the  nighttime 
topside  ionosphere  at  low  and  moderate  solar  activity,  it  was  discovered  that  a 
meridional  neutral  wind  could  produce  a  trough  in  the  ion  temperatures  that  re¬ 
sults  in  the  ion  temperature  being  cooler  than  the  expected  neutral  temperature 
over  portions  of  the  field  lines.  This  trough  in  the  ion  temperature  was  found  not 
to  form  in  the  altitude  region  where  there  is  strong  interaction  between  the  ions 
and  neutrals  because  the  neutrals  act  to  quench  any  heating  or  cooling  effects 
exerted  on  the  ions.  The  latitudinal  extent  of  this  trough  is  determined  by  field 
line  geometry  and  the  characteristics  of  the  F-region  neutral  winds  [Bailey  and 
Heelis,  1980]. 

Interhemispheric  flow  has  been  shown  to  significantly  affect  the  plasma  tem¬ 
peratures,  and  that  for  solstice  conditions  the  flow  of  pleisma  is  generally  from 
the  sununer  to  the  winter  hemisphere.  This  interhemispheric  flow  results  in 
warmer  temperatures  in  the  winter  hemisphere  due  to  adiabatic  compression  as 
the  plasma  descends  along  the  field  line  [Bailey  et  al.,  1975;  Bailey  et  al.,  1973]. 
The  electrodynamic  drift  has  little  effect  on  plasma  temperatures  at  night  [Bai¬ 
ley,  1983;  Bailey  et  al.,  1973],  but  is  significant  during  the  day  [Bailey,  1983]. 
Using  the  coupled,  time-dependent  model  of  Bailey  [1983],  Moffett  et  al.  [1986] 
demonstrated  that  for  moderate  solar  activity  and  solstice  conditions  thermal 
coupling  produces  very  strong  upward  flows  of  O"*"  in  the  topside  summer  iono¬ 
sphere  at  the  time  of  sunrise  in  the  winter  hemisphere  and  a  lesser,  but  opposite, 
effect  at  the  time  of  winter  sunset.  For  these  conditions,  counterstreaming  be¬ 
tween  the  and  O'*"  ions  occurs  with  the  total  flux  at  500  km  being  primarily 


composed  of  !!■*■  ions.  Additionally,  the  ions  were  found  to  exhibit  strong 
upward  and  downward  flows  at  sunrise  and  sunset,  respectively. 

Since  the  ions  are  primarily  heated  by  the  electrons,  it  is  appropriate  to  now 
discuss  some  of  the  properties  of  the  electron  gas.  Extreme  ultraviolet  radiation 
is  the  primary  energy  source  for  the  electrons  during  the  day.  Energetic  pho¬ 
toelectrons  are  created  through  the  photoionization  of  the  neutrals  and  these 
photoelectrons  lose  their  energy  in  elastic  collisions  with  the  ambient  electron 
gas  and  inelastic  collisions  with  the  neutrals  and  ions  [Schunk  and  Nagy,  1980]. 
At  the  heights  being  considered  by  our  model,  the  primary  loss  to  the  neutrals 
occurs  through  the  fine  structure  excitation  of  O  [Bailey  and  Heelis,  1980).  Be¬ 
low  250  -  300  km,  photoelectrons  lose  their  energy  locally;  above  this  height, 
they  escape  along  magnetic  field  lines,  heating  the  ambient  electron  gas  at  high 
altitudes.  This  results  in  the  electron  temperatures  being  greater  than  the  ion 
temperatures  during  the  day  [Schunk  and  Nagy,  1980). 

Model  and  Assumptions 

This  thesis  presents  the  results  of  a  time-dependent,  self-consistent  model  of 
the  equatorial  ionosphere.  This  model  is  intended  to  be  incorporated  with  other 
models  covering  other  ionospheric  regions  (mid  and  high  latitudes)  to  form  a 
three-dimensional,  global  model  of  the  terrestrial  ionosphere.  In  this  model, 
the  continuity,  momentum,  and  energy  equations  for  O"^,  H"^,  and  electrons 
axe  coupled  and  solved  in  a  quasi-simultaneous  manner  by  first  solving  for  the 
coupled  H"*"  and  densities  and  then  solving  for  the  electron  temperature,  the 
H'*’  temperature,  and  the  O'*"  temperature  at  each  time  step.  For  the  daytime 
case,  we  use  an  empirical  expression  for  the  heating  rates  due  to  photoelectrons 
which  were  obtained  from  the  results  of  Swartz  et  al.  [1975]  (to  be  discussed  in 
more  depth  in  Chapter  III).  Our  neutral  atmosphere  comes  from  the  MSIS86 
empirical  model  [Hedin,  1987). 
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In  the  formulation  of  this  equatorial  model,  a  number  of  assumptions  were 
made  which  will  now  be  commented  on.  The  electron  densities  are  calculated 
through  the  assumption  of  charge  neutrality,  and  the  field-aligned  velocities  for 
the  electrons  are  obtained  by  assuming  that  there  are  no  currents  flowing  parallel 
to  B.  This  assumption  is  supported  by  the  results  of  Kendall  and  Pickering 
[1967]  who  found  that  for  reasonable  estimates  of  the  parallel  current  in  the 
low-latitude  ionosphere,  the  effects  of  this  parallel  current  c'\  the  field-aligned 
drift  velocity  is  negligible  and  that  the  standard  view  of  ambipolar  diffusion 
along  field  lines  is  correct.  We  consider  an  ionosphere  composed  of  two  ion 
species,  H'*'  aind  O'*',  and  electrons.  This  assumption  is  justified  because 
and  H'*'  are  the  dominant  ions  between  the  F2  peak  and  3000  km,  and  NO"^ 
ceases  to  be  the  dominant  ion  between  150  and  250  km  [Banks  and  Kockarts, 
1973;  Anderson  and  Rusch,  1980).  Additionally,  Sehunk  and  Walker  [1973]  have 
shown  that  the  density  is  not  significantly  affected  by  £^-region  polaurization 
electric  fields  and  molecular  ion  densities  in  the  region  where  is  the  dominant 
ion.  Our  use  of  the  fluid  equations  assumes  the  ionospheric  gas  mixture  has  a 
Maxwelliam  velocity  distribution  and  that  along  the  field  lines  the  mean-free-path 
for  collbions  between  the  H"*"  (neutrals  at  low  altitudes)  and  O'*"  ions  is  smaller 
than  the  scale  height;  these  are  true  in  the  equatorial  ionosphere.  The 
fluid  behavior  of  the  plasma  is  also  enhanced  by  the  presence  of  the  magnetic 
field  which  acts  to  create  plasma  flows  parallel  to  B  [Chen,  1984).  For  diffusion 
parallel  to  B  to  exist,  the  gyrofrequency  of  a  species  must  be  greater  than  the 
collision  frequency  for  the  species  [Banks  and  Kockarts,  1973).  This  condition 
is  met  at  all  altitudes  being  examined  by  our  model,  except  at  the  very  base  of 
the  field  line  where  photochemistry  is  the  dominant  process. 

For  the  last  10  years,  a  number  of  models  of  the  low-latitude  ionosphere  have 
taken  into  account  the  electrodynamic  drift  of  the  plasma.  However,  these  stud¬ 
ies  involved  very  “idealized”  electric  fields  or  assumed  a  value  for  the  equatorial 
drift  velocity  [Bailey,  1983;  Bailey  and  Moffett,  1979;  Sterling  et  al.,  1969).  In 


some  cases,  the  electric  field  or  the  electrodynamic  drift  velocity  was  chosen 
to  obtain  the  best  fit  to  a  set  of  observations  [Anderson,  1973o,  19736].  Only 
Chan  and  Walker  [1984a,  19846]  used  an  electric  field  model  that  was  based 
on  observations  over  a  variety  of  conditions  at  more  than  one  observing  site. 
Unfortunately,  they  assumed  the  vertical  drift  velocities  were  invariant  along  a 
given  field  line.  This  is  probably  an  erroneous  assumption,  since  field  lines  are 
generally  considered  to  be  equipotential  [Richmond  et  al.,  1980]  and  dipole  field 
lines  converge  as  one  moves  downward  in  altitude  along  them.  This  should  lead 
to  a  variation  of  the  electric  field  along  the  field  line,  and,  hence,  a  variation 
in  the  electrodynamic  drift  velocity.  We,  therefore,  calculated  the  variation  of 
electrodynamic  drift  velocity  along  a  given  field  line  by  exploiting  the  equipo- 
tentiality  of  the  field  line  (See  Appendix  C  for  details  of  our  method).  Using 
these  techniques,  we  feel  the  variation  of  the  electrodynamic  drifts  along  the 
field  line,  with  season,  and  with  local  time,  avoid  a  priori  assumptions  concern¬ 
ing  densities  along  the  field  line  or  electric  field  variations.  This  should  increase 
the  general  accuracy  and  applicability  of  our  model. 

Additionally,  few  studies  [Anderson,  1973a,  19736;  Chan  and  Walker,  1980a, 
19806;  Young  et  al.,  1979]  have  treated  the  dipole  field  as  having  a  tilt  in  re¬ 
lation  to  the  earth’s  rotational  axis,  except  perhaps  through  incorporating  the 
effects  in  the  “idealized”  neutral  winds  that  were  assumed.  Our  model  provides 
self-consistent  values  for  the  neutral  winds  as  a  function  of  latitude,  altitude, 
season,  solar  activity  (maximum  or  minimum),  longitude,  auid  local  time.  De¬ 
tails  of  our  neutral  wind  model  are  included  in  Appendix  A.  It  should  be  noted, 
however,  that  our  neutral  wind  model  is  empirical  and  the  coupling  of  both  an 
ionospheric  model  and  a  neutral  atmosphere  model  is  required  for  truly  consis¬ 
tent  calculations.  It  is  felt  that,  like  the  electric  field  model,  this  neutral  wind 
model  enhzuices  the  general  applicability  of  our  overall  modelling  efforts. 

Finally,  it  should  be  noted  that  while  the  azimuthal  component  of  the  elec¬ 
trodynamic  drift  hzis  been  assumed  negligible  based  on  the  results  of  Sterling  et 


CHAPTER  II 


CONTINUITY  AND  MOMENTUM  EQUATIONS 


An  Orthogonal,  Curvilinear  Coordinate  System 


We  wish  to  model  the  variations  of  ion  and  electron  temperature  and  density 
for  a  topside,  low-latitude  ionosphere  composed  of  two  ion  species,  O"*"  and  H"^. 
To  accomplish  this,  the  plasma  is  assumed  to  diffuse  parallel  to  the  magnetic 


field  line  and  drift  normal  to  the  field  line  due  to  E  x  B  forces.  For  a  magnetic 


dipole,  it  is  convenient  to  adopt  a  curvilinear,  orthogonal  coordinate  system 
which  has  one  coordinate  parallel  to  B  and  one  coordinate  normal  to  B.  The 
coordinates  in  this  system  may  be  defined  as  p,  q,  and  <j>  where  ^  is  the  magnetic 
longitude.  The  location  of  any  point  along  a  given  field  line  may  be  uniquely 
defined  in  terms  of  p  and  q  where  the  field  lines  are  represented  by  a  family  of 
curves  along  which  p  is  constant,  and  the  equipotential  surfaces  are  represented 
by  the  family  of  curves  along  which  q  is  constant.  Expressions  relating  the  p-q 
coordinate  system  to  spherical  coordinates  are 


Te  sin'  9 


r|  cos  9 


where  Vg  is  the  radius  of  the  earth,  r  is  the  radius  of  any  given  point,  and  9  is 


the  magnetic  colatitude.  The  equation  for  a  field  line  is 


r  =  a  sin^  9 


where  a  is  the  radius  of  the  field  line  at  the  equator.  Since  we  are  interested  in 


parallel  motion,  the  unit  tangent  to  the  field  line  is  defined  as 


i  =  (sin  /,  cos  /,  (f>) 
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with  i  being  defined  positive  above  the  horizontal  and  directed  towards  the 
geomagnetic  equator  in  the  northern  hemisphere.  If  s  is  a  unit  arc  length  along 
the  field  line,  then 

d  .  d  cos  I  d  d 

ds  dr  r  dO  dq 

where  tr  =  +  3cos^  [Kendall,  1962].  See  Figure  1  for  an  example  of  Jiow  q 

varies  along  a  field  line. 


Colatitude 

Figure  5.  A  diagram  showing  the  variation  of  q  along  a  magnetic  field  line  that 
crosses  the  equator  at  3000  km  altitude.  The  solid  curve  represents  the  field  line 
(or  a  curve  of  constant  p)  and  the  dashed  curves  show  lines  of  constant  q.  Aq  is 
chosen  to  be  0.1. 


Transport  Equations 


At  low  and  middle  latitudes  the  ion  or  electron  number  density  is  determined 
by  two  main  processes.  These  are  atmospheric  photochemistry  and  transport. 
Both  processes  cam  be  effectively  modelled  using  the  equations  for  the  conser¬ 
vation  of  mass  and  momentum  so  long  as  the  plasma  can  be  thought  to  behave 
as  a  fluid.  These  conservation  equations  are  obtained  by  taking  the  moments  of 
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the  Boltzmann  equation.  The  continuity  equation  can  be  written  as 


dn. 


+  div  (n«Ua)  =  Pi  -  LiUi 


(6) 


where  n,  is  the  number  density,  Ug  is  the  plasma  drift  velocity  for  the  species, 
Pg  is  the  photochemical  production,  and  Lg  is  the  chemical  loss  rate  The 

momentum  equation  is  of  the  form 


ngnii 


+  VP*  +  S7  -U-  ngmgG 


-  nge. 


Ug  X  B 


^Mg 

St 


(7) 


where  m,  is  the  mass  of  species  s,  p,  =  rigkTg  is  the  partial  pressure  for  an  ideal 

D  /I 

gas,  Cg  is  the  charge,  Tg  is  the  temperature,  Tg  is  the  stress  tensor,  TTt  - 
is  the  convective  derivative,  G  is  the  acceleration  due  to  gravity,  E  is  the  electric 
field,  B  is  the  magnetic  field,  c  is  the  speed  of  light.  Mg  is  the  rate  of  momentum 
exchange  due  to  collisions  of  s  with  other  species,  and  k  is  Boltzmann’s  constant 
[Schunk,  1983]. 


Atmospheric  Chemistry 


Although  this  model  is  primarily  valid  for  the  diffusion  dominated  portion  of 
the  ionosphere  where  charge  exchange  between  and  O'^  are  the  dominant 
chemical  reactions,  the  ion-atom  interchange  reaction.^  of  O"*"  with  O2  and  N2, 
and  the  photoionization  of  atomic  oxygen  are  also  included. 

Chemical  reactions. 

0+ 

0+  +  NO-^  +  0,  ki  (8) 


ki  =  1.533  X  10“*^  -  5.92  x  10"*^—  -t-  8.6  X  10 


\300j 


300°  <T  <  1700°  K 


{9a) 
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ky  =  2.73  X  10“*^  -  1.155  X  10"‘* - f-  1.483  X  10~‘® 


300 


\300j 


1700°  <T  <  6000°  K 


(96) 


0+  +  O2-*  0+  +  0,  k2 


(10) 


ki  =  2.82  X  10“”  -  7.74  X  10“^*—  +  1.073  x  10"^^  ( —  J 

300  V300 / 


-5.17x10“^^/^^)  +9.65  X  10"'®  f—) 

VSOO/  V300/ 

300°  <  T  <  6000°  if 


(11) 


H+ 


0+  +  H  -*  H+  +  0,  k3  =  2.5  X  10"" 


Production  :  O'*'  +  if  — >  H'*'  +  O,  kz 


(12) 


(13) 


Loss  :  H+  +  O  -*  H  +  0+  ,k4  =  2.2  X  10““rj^''+  (14) 

The  reaction  rates  are  in  cm^s"'  and  the  temperatures  are  in  °K  [5cAunA:,  1983]. 

Photoiomzation.  At  low  altitudes,  the  primary  source  for  is  through  the 
photoionization  of  neutral  O.  The  Chapman  production  function  describing  this 
reaction  is  of  the  form 


P{.Kx)  =  P<x>n{0)ex^[-a[h,x))  (15) 

where  Poo  is  the  total  production  of  ion-electron  pairs  for  an  overhead  sun  over 
all  ionizing  wavelengths  of  O,  h  is  the  altitude  of  the  point,  x  is  the  solar  zenith 
angle,  and  a  is  the  optical  depth.  Poo  is  based  on  the  integrated  solar  flux  values 
in  seven  bands  from  165A  to  91lA  and  the  average  ionization  cross-sections  of 
O  for  each  band  assuming  each  absorbed  photon  produces  an  ion-electron  pair. 


The  value  of  1.1  X  10  ^5  ^  was  adopted  for  solar  minimum  conditions  [Banib 
and  Kockarts,  1973).  The  variation  of  jPoo  with  sunspot  number  is  given  by 

Poo  (P)  =  1-1  X  10“^  (1  +  0.0124P)  (16) 


where  R  is  the  average  sunspot  number.  The  optical  depth  is  given  by 

«(fc,x)  =<’AO)n(0)H(0)Ck(^j~,x'^  +<T.(JVj) 

•n(Wi)/f(JV2)C/.(^^,x)  (17) 

where  Oa  is  the  ionization  cross-section,  H  is  the  scale  height,  and  Ch  (P’X) 
is  the  Chapman  function.  The  adopted  values  for  the  ionization  cross-section 
are  <Ta  (O)  =  0.9  x  10~^^  and  {^2)  =  1-5  x  10~*^  [Anderson,  19736).  The 
equations  for  the  Chapman  function  for  any  point  p  take  into  account  grazing 
incidence  and  are  as  follows: 

For  Xp  <  90° 


c  +  dy  +  y^' 

^  f 


0  <  y  <  8 


8  <  y  <  100 


(18) 

(19) 


For  Xp  >  90° 


Ch  {Xp,  Xp)  =  (27rXp)5  [sin  (xp)^  exp  {Xp  (1  -  sin  Xp)) 
1  a  -h  6y  1 


(20) 


Ch  {Xp,  Xp)  =  (2jrXp)^  [sin  (Xp)^  exp  (Xp  (1  -  sin  Xp)) 
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¥ 

f 

4 


where 


y=(^)  Icosxpl, 


~  'Jj' 

a  =  1.0606963,  b  =  0.55643831,  c  =  1.0619896,  d  =  1.7245609,  /  =  0.56498823, 
and  g  =  0.06651874  [Smith  and  Smith,  1972].  The  solar  zenith  angle  was  com¬ 
puted  using  the  formula  found  in  Schunk  and  Walker,  [1973]  which  is 


cos  Xp  ~  sin  ^  sin  <p  cos  ^  cos  ip  cos  ^ 


27r  {HL- 12)' 


where  ^  is  the  solar  declination  angle,  HL  is  the  local  solar  time  in  hours,  and 
(p  is  the  geographic  latitude.  From  the  above  equations,  the  dominamt  photo¬ 
chemical  reactions  can  be  accounted  for  with  the  exception  of  the  reactions 
involving  NO"*"  which  are  only  important  at  the  lowest  altitudes  being  considered. 

Momentum  Equation 

The  calculation  of  the  time-varying  ion  and  electron  number  densities  along 
the  field  line  requires  the  parallel  component  of  the  ion/electron  drift  to  be 
specified.  These  parallel  drift  velocities  can  be  obtained  from  equation  (7).  First, 
if  the  following  assumptions  are  made:  (a)  steady  state  =  O),  (b)  the  parallel 
drifts  are  subsonic,  therefore,  the  nonlinear  terms  in  (7)  can  be  neglected,  (c) 
heat  flow  in  the  collision  term  is  negligible,  (d)  no  field-aligned  currents, 

(e)  flow  is  parallel  to  B  (u,  X  B  =  0),  and  (f)  viscous  effects  can  be  neglected, 
equation  (7)  becomes 


VP<  -  n,mtG  -  n,etE  = 


Now,  by  assuming  a  perfect  gas  and  realizing  that  along  the  field  line  G  —  g  sin  It 
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where  i  is  a  unit  vector  parallel  to  B,  equation  (24)  can  be  rewritten  as 


—  {ntkTt)  t  +  sin /<  —  n«ejE  =  — 

OS  dt 


(25) 


The  electric  field  in  equation  (25)  can  be  obtained  from  the  electron  mo¬ 
mentum  equation  by  making  the  same  assumptions  as  before  and  by  neglecting 
gravity  and  the  collisional  terms  due  to  the  much  smaller  mass  of  the  electrons. 
This  gives 


E 


(26) 


Substituting  equation  (26)  into  (25)  and  assuming  charge  neutrality  gives 


kTt%~t  +  risk 


dT,, 


ds 


-  t  +  ntTn,g  sm  It  +  krig-^t 
as  as 


^  kTeUgdNg^  _  ^M,  _  ^ 


Ng  ds 


6t 


(27) 


Next,  the  rate  of  transfer  of  momentum  due  to  collisions  needs  to  be  evalu¬ 
ated,  Since  general  expressions  for  this  quantity  are  only  available  for  limited 
particle  distributions,  it  becomes  necessary  to  assume  the  ionosphere  approxi¬ 
mates  a  displaced  Maxwellian  distribution.  This  leads  to  the  5-moment  approach 
since  it  involves  density,  temperature,  and  three  components  of  velocity  for  each 
species.  Ignoring  the  heat  flow  ternos,  as  has  already  been  assumed,  the  rate  of 
momentum  transfer  can  be  expressed  as 

~  ^  ("s'  -  "t)  (28) 

where  Ugt  is  the  collision  frequency  for  momentum  transfer  between  species  s 
and  t  (5chunA:,  1983|. 

Substituting  equation  (28)  into  (27)  and  rewriting  the  resulting  equation  gives 

„  t/gt _ kTg  /  1  drig  IdTg  m,gsin/ 

^*'i'gn +  t^st  mg{ugn  + i'gt)\ng  ds  Tg  ds  kTg 
1  dTg  Tg  dNg  VgnU  cos  /m, 

¥,~^^  TgNg^s  fcn 


(29) 


where  U  cos  It  is  the  parallel  (to  B)  component  of  the  meridional  neutral  winds. 


If  we  define 


then  the  diffusion  coefficient  becomes 


D,  =  — ^A.. 


By  also  defining, 


Q  _  ( '^*9  sin  /  1  dTg  1  dT^  cos  I\  « 


the  momentum  equation  becomes 

1/  « *  r,  /  1  dn,  Tt  dNc  ^  \ 

U"  -  +  3^17 

Equation  (33)  is  now  ready  to  be  substituted  into  the  continuity  equation. 
Continuity  Equation 

The  continuity  equation  allows  us  to  calculate  the  time- varying  ion/electron 
densities  along  the  field  line.  Using  Sterling  et  al.  ’s  [1969]  technique,  Ug  can  be 
broken  down  into  its  components  giving 

u,  =  Uem  +  Uth  +  Ug  (34) 

where  uth  is  the  vertical  velocity  of  the  neutral  air  due  to  thermal  expansion  and 
contraction,  and  Uen,  is  the  electrodynamic  drift  of  the  ionosphere.  Equation 
(6)  now  becomes 


=  Ps~  L,ng  -  div  (n,u<,n,)  -  div  (n,Ug)  -  div  (n*U(fcsin  It)  .  (35) 
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Now,  using  the  following  definitions: 


-div  {rifUtij)  =  i  ■  grad {—ngUth)  —  UgUthdivt 


(36) 


[Sterling  et  aK,  1969], 


div  (t)  = 


9  cos  0  +  15  cos^  0 
^  ’ 


(37) 


i  -  grad  = 


ds 


(38) 


[Kendall,  1962],  and  equation  (33),  the  various  divergences  can  be  expanded. 

In  general,  the  vertical  velocity  of  a  surf2u:e  of  constant  pressure  cam  be  written 


as 


n  dToo  {z  —  Zo)  r 
Too  d<j>  Te  +  Zo 


(39) 


assuming  T  =  Too  and  Too  is  independent  of  z  where  Zo  =  120  km  atnd  fl  is  the 
earth’s  rotation.  Now,  using  equations  (36)  -  (39) 


div  (— n^Ufftsin/t)  = 


2  cos  Or^uth  dn. 


+  ngAz 


d  {uth  sin  /) 


—  ngUth  sm 


in  /  ^ 


dq  ‘  ‘  dq 
9  cos  0  +  15  cos^  0 


ra 


zjz 


) 


(40) 


The  electrodynamic  drift  can  also  be  separated  into  its  components  giving 


u«ni  =  Uj.  +  +  nr  sin  0^. 


(41) 


Since  the  div 


— n^flrsin  0 


=  0, 


div  (-n,Ueni)  =  div  (-n,uj.)  +  div  (-n,u^)  . 


(42) 


Recall,  that  we  are  interested  in  the  drifts  and  densities  along  a  given  field  line, 
and  that  it  is  convenient  to  allow  the  coordinate  system  to  move  with  velocity. 
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Uem-  Using  the  definition  for  the  convective  derivative,  we  get 


Dng  _  dn, 
Dt  ”  ~dr 


The  divergences  of  the  electrodynamic  drifts  are 


div  (u j.)  =  ■  ^  (6  cos®  0  —  Z  cos*  0  —  A  cos^  0  +  l)  (44) 

arc  /  V  / 


where  uj.o  is  the  value  of  Ux  at  the  equator. 

Substituting  equation  (33)  into  (35)  and  using  previous  assumptions  give 


j.-..  /  _  ^  ^2  dnc  ,  /  Al  dnt  dn,  ,  A\  d^ric 

d,v  (-n.u.)  -  —g~-—  + 

-  MAzdn,\  dWg  rtgAz  dNc  ^  rAtUg  d^Nc 

nt  dq  )  dq  ntNg  dq  *\ntNc  dq^ 
Az  dNc  dn,  Azrig  dNc  A^ng  / dNc  \ 

^  NcTH  dq  dq  Ncn\  dq  dq  ntN^  \  / 

_  Azrig  dNc\  dVg  A^rig  _y  fAz  drig 
ntNe  dq  /  dq  nt  "\ntdq 
Azn-g  dnt  ^  ng 


where  nt  is  the  number  density  of  the  other  ion,  L/g  =  Dgnt,  Vg  = 


Az  =  — 


fVg  = 


AzTcD' 


•4.-1  — 


9cos0  +  15  cos^  0 


Now,  by  combining  equations  (35),  (40),  (42),  and  (44)  -  (46)  we  can  write  the 
general  form  of  the  nonlinear,  second-order  partial  differential  equation  which 
couples  both  ions  and  gives  the  time-varying  ion  densities  along  a  field  line  as 


=  Pt  —  Lftit  —  fig  ^6  cos®  0  —  Z  cos'*  0  —  A  cos*  0  -f  j 

1  2cos0r^utkdng  3(u/*sin/) 


dWg  A2ng  dN, 


A^Az  an 


d 

A^ng  dnt 


ntNg  dq*  Ngnt  dq  dq  NgU^  dq  dq 


Azng 

ntm 


ig  fdNcV  A3ngdNg\ 

7  \  ^0  J  ntNc  dq  )' 


To  solve  this  equation  numerically,  it  becomes  necessary  to  linearize  all  terms. 
This  is  accomplished  by  performing  a  time-varying  Taylor’s  expansion  on  the 
nonlinear  terms  and  then  substituting  them  back  into  equation  (50).  After 
linearization,  a  second-order  differential  equation  for  the  ions  of  the  form 

^  ^  dug  ^  ^  d^nt  ^  dnt  ^  . 

results.  The  various  coefficients  for  each  species  are  presented  below. 


6*2  =  A2u"ff^A2  +  A25-  {A2D2)  (  — - —  A2AzD'2  ( — - ^ 


Ca  =  2.2  X  (O)  -  (A.D',)  j 


I  ^  dNQ+  dNfj+\  2p» 


1 

d^No.\ 

N^ 

dq^  J 

+ 


+ 


.  A  r.>(  1  dNoA°  A  9^2  f  No^  dN. 


Nq^  dN^Y  ,A^( ^ 
NeN^.  dq  )  ^  39  I  Nl+ 


+  -^2^2 


dNp+y 


—  A2V2 


2No.aN„,y 


^H+  ) 


'•0+ 


H*  J 

No+  d^N, 


^H+ 


No.  a^N^Y  /  1  dNo^dN. 

Nu.N?  dq^  J  ^  dq  dq 


+ 


Nh^N!  dq^ 

1  dNo+  dN 


+ 


^H+^e  dq  dq 

2No+  dNtf+dNc 
NeN^+  dq  dq 

2No.  (dNc^^ 
Nh.N!  V  dq 


)••«■•(!£»  O’ 


+ 


N^-dq^J  dq 

^  A  ^*^2  /  Nq.  dN,y  .  ^  fNo.Y  _  .  „  fj_d^\ 
*  39  VN^+Ne  3?  y  *  3g  \Nff+J  ^  *  \-^H+  ^9  / 

,  * «...  (feO' 


59  J 


>l2l^2 


JVqa.  dNfj.dN, 
^NcN]j+  dq  dq 


(  No.  fd^yX 

{Nfj.N^  \dq  }  ) 


+  f’oon(0)exp(-a(/i,x)) 


(58) 


.  ;  " 


0.077Nh+ 
t'o+H+  -  -rsTi 

^0+H+ 


=  6.61  X  10-“n  {H)  Till  (1  -  0.047 logic  Th+)"  +  1-32  X  lO'^^n  (He) 
+  3.67  X  10-^^n  (O)  (l  -  0.064  logic  Trf 
+  6.82  X  10-‘°n  (Ni)  -t-  6.64  x  lO'^^n  (O2) 


ft  =  Aju"  *  Ai  +  A,^  (A,D;)  (^)  -  AjAjDl  [-j^) 


~AzWi 


C3  =  -2.2  X  10-“4?n(O)  -  ^ 

\  /  1 

_  (6  cos®  0  -  3  cos'*  0  -  4  cos*  0  -f  1 )  -  ) 

d  N  atVi  /  1  ^  _  JiHi_^y 

(’*0+^0  +  ^*"^VNo+JVc  59  No+^l  ^9  ^ 

. ...  /  1 

+  ^*^*VHo+N«  a?*  Ho+^^e*  ^9*  ^ 

/  .1  dN„+dN,y  ^ 

a,  a? 
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Nfj+  dNQ+  dN, 


N^+N^  dq  dq 
2Nh+ 


_  (^YY  _  .  ^  (  1  dN, 

^i[dq))  dq 


No^Ni 
Njf^  aN,y 


No+N^  dq 

“0+^l-^3 


\  °  d 

j  +>42^(utAiSm/)  -  ylsatAsin/ 


(63) 

(64) 


N; 


H+ 


AiWi 

A^Wy 


(  1  dNji^Y  ,  /jv 

VNo+N,  dq  J  ^  '[N, 

(  2Nh+  dN^Y 
\No^N^  dq  J 


Nh^  dNo^  ^  Njf+  dN, 


A3W1 


Nl^  dq  N,Nl^  dq 

(J^Y 


\No^Nj 


(65) 


2.5 . 10-“ryv  (H)  -  A,±  (MB’,) 


2  BN„taNo*y  1 


dq  dq 


-AlV, 


Nl+  dq^ 

dN, 
>>  dq 


+ 


Nh^ 

dN,\ 

dq  / 

A2Vy(^ 

2Nff+ 

Nh^ 

B^N, 

+  ^42^ 


1  dNff^ 


dq 


-  A^Vy  ^ 

NhA 

0 

-A^Wy 

(  1 

+ 

dN, 

dq 

dq 

dq  J 

Nh+  d^N, 


N,N^^  dq^ 


1  dNff.dN,Y  /  dNo^dN 

Nq+N}  dq  dq  ) 

2N„*  dNpt  dN,y  (  N„,  /a^y 

N,Ny  dq  Bq)  ' /V*+ /Vf  (  3?  j 


% 

s 

> 


29 


No^Nf  \dq  J  )  ^  '\No+N^  dq  N,Nl^  dq  ) 


(66) 


C7  =  A2 


dq 


dq  j 


iVo+  dq^ 
dy^(_Nfj+  dN, 
dq 


No+N,  dq 


'  )  dq  \NoJ  '  'Uo+  dq  ) 


+ 


+ 


dq  dq 


Njj+  dNp+dN^' 

cN^+  dq  dq 


I  ffH* 


-.(ifef)' 


(67) 


^H+0+  = 


I.23N0+ 

m3/2 

^0+H+ 


(68) 


^H+n  =  2-65  X  10"^°n  (H)  (1  -  0.0831  log,o  TV)^  +  10.6  x  lO'^^n  (He) 

+  6.61  X  10““ n  (O)  (1  -  0.047  logjo  7^+)^ 

+  33.6  X  10“^®n  {N2)  +  32.0  x  10“^°n  (O2)  (69) 


Electrons 

Assuming  charge  neutral'  ,  ,  the  electron  number  density  is  given  by 


N,  =  No^+N„+.  (70) 

In  addition,  the  assumption  that  there  are  no  currents  flowing  parallel  to  B 
allows  the  parallel  drift  velocity  for  the  electrons  to  be  calculated  from  the  drift 
velocities  and  number  densities  of  the  ions.  The  expression  for  the  parallel, 


4^  *»•  4^  iy> J^»  !■*  ^*'N* 


electron  drift  velocity  is 


Ne 


^  iiAb  k_-\«  ,  -\.  ,  K.  .  *  .^  - 
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CHAPTER  III 


ENERGY  EQUATIONS  FOR  0+,  H+,  AND  ELECTRONS 


The  plasma  heat  balance  equation  for  O"^,  H'*",  and  electrons  can  be  derived 
from  the  general  form  of  the  energy  equation, 

Wt  (^)  (V  •  Us)  +  V  •  Qb  +  rg  :  VUb  =  ^  +  Qs  -  Lg  (72) 

where  ^  ^  +  Ug  •  V»  the  convective  derivative  of  species,  s;  p,  =  rigkT,, 

the  partial  pressure  of  species,  s;  rig  is  the  number  density;  k  is  Boltzmann’s 
constant;  Tg  is  the  temperature;  Ug  is  the  drift  velocity;  Qg  is  the  heat  flow 
vector;  Tg  is  the  stress  tensor;  Qg  is  the  heating  rate;  Lg  is  the  cooling  rate; 
V  is  the  coordinate-space  gradient;  and  is  the  rate  of  energy  exchange  due 
to  collisions.  As  in  the  continuity-momentum  equation,  we  assume  a  displaced 
Maxwellian  for  the  plasma.  The  approximate  expression  for  the  rate  of  energy 
exchange  due  to  collisions  for  the  5-moment  approax:h  takes  the  form 

^Eg  3/:  (Tj  —  Tg)  -f  m(  (ug  —  Ug)^  $,(  ^  ^ 

-rr  =  ngmg  y  Ugt - — — -  17.JJ 

ot  "  rrig  +  mt 

where  ’J^gf  and  $,4  are  velocity  dependent,  correction  factors  and  Ugt  is  the 
collision  frequency  for  momentum  transfer  for  species  s  and  t. 

Neglecting  viscous  effects  and  using  the  continuity  equation  in  the  form 

^  =  -div  (rigUg) ,  (74) 

which  is  valid  in  a  diffusion  dominated  region,  and  the  definition  of  qg, 

Qb  —  “-^B  V  Tg,  ('^^) 

where  Xg  is  the  thermal  conductivity  [Schunk,  1983],  allows  equation  (72)  to  be 
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written  in  the  form 

3  DT 

=  Q'»  -  A:n,r,<ftt;ug  +  div  (Ag  V  Tg)  +  Fg,o  (76) 

where  Q'g  is  the  net  rate  of  heat  gain  due  to  collisions  and  Ff^o  is  the  frictional 
heating  due  to  relative  field-aligned  motion  between  an  ion  and  the  neutral  gas 
[Bailey  et  al.,  1973].  The  drift  for  a  given  species  is  again  divided  into  its  parallel 
and  perpendicular  components  such  that 

Ug  =  u"  -H  Uem  (77) 

where  u"  is  the  parallel  drift  and  is  the  electro  dynamic  drift  [Sterling  et  al., 
1969).  This  frictional  heating  term  has  the  form 

Fs,o  =  n,Vt,o  (wg  -  U  cos  /)^  (78) 

where  i/g,o  is  the  collision  frequency  for  momentum  transfer  between  the  ion  and 
neutral  species,  U  is  the  meridional  component  of  the  neutral  wind  velocity,  and 
I  is  the  magnetic  dip  angle.  The  collision  term  for  O'*"  is 

I'O+O  =  3.42  X  10““n  (O)  (1.04  -  0.064  logjo  Tr)^  [s“^]  (79) 

and  for  it  is 

^H+O  =  6-61  X  10'^‘n  (O)  (1  -  0.047  log^o  Trf  [s"^]  (80) 

[Bailey  and  Heelis,  1980]  where  Tr  {To+  +  Tn)  [Raitt  et  al.,  1975].  Since  the 
primary  conduction  of  heat  is  parallel  to  B, 


where  B  is  the  earth’s  magnetic  induction  emd  s  is  the  distance  along  the  field 
line  [Bailey  et  al.,  1973].  Now, 


div  (Ug)  =  B^  +  div  (u* 


-v 


-  ■urj’ 


and  div  (u*m, )  is  given  by  equations  (44)  and  (45) .  We  again  assume  the  co¬ 
ordinate  system  is  drifting  at  Uem-  Now,  the  time-dependent  form  of  equation 
(76)  b 

•  ^6  cos®  6  —  3  cos^  6  —  4  cos^  6  +  1^ 


rsinp  op  J  as  B  as 

[Bailey  et  al.,  1973;  Bailey  and  Heelis,  1980].  Using  the  definitions  for  B, 

=  (84) 

where  Bo  is  0.32  xl0~^  Teslas  and  R  =  ^  [Hargreaves,  1979], 

d  .  cosi  d  .  . 

^  =  +  — a«’ 

[Baxter,  1964],  and  by  realizing  variations  in  the  temperature  of  a  given  species 
in  the  r-direction  are  incorporated  into  the  variations  in  the  ^-direction,  equation 
(80)  can  be  expressed  as 
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where  A,  =  [Bailey  et  al.,  1973].  Next,  using  the  definition  of  a  dipole 

field  line,  r  —  r^  sin*  6  and  previous  assumptions,  the  field-aligned  component  of 
the  second  term  on  the  right-hand  side  of  equation  (83)  becomes 


•  rr, 

-knsTsB——  = - r- 

ds  B  TeSin 


.r,  /a< 

V  96 


3  sin  26u'l 


+  6  cot  6u. 


Substituting  (84)  -  (87)  into  equation  (83),  produces  a  nonlinear,  second-order 
differential  equation  of  the  form 


where 
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and  for  ions 


or  for  the  electrons 

Frictional  heating  due  to  collisions  with  the  neutrals  is  ignored  for  the  electrons. 
The  expressions  for  the  heat  transfer  terms  are  given  below.  For  the  ions,  this 


expression  is 


Q'o+  -  Qo+e  +  Qo+,H+  +  Qo+n  +  *5o+,if 


Q'h+  -  Qn+e  +  Qh+.0+  +  Qla+n  + 

where  Qae  is  the  heating  rate  due  to  elastic  collisions  of  ion  species,  s,  with  the 
electrons,  Qat  is  the  heating  rate  due  to  elastic  collisions  between  ion  species  s 
and  t,  QjJj  is  the  heating  rate  due  to  elastic  collisions  with  the  neutrals,  and 


Q*”  is  the  heating  rate  due  to  inelastic  collisions  with  the  neutral  oxygen  or 
hydrogen.  For  the  electrons,  is 

Qe  —  QeO+  +  QeH+  +  Qf.st.O  "t-  Qh  (96) 

where  Qeo+  electron  heating  rate  due  to  collisions  with  O"^,  QgH+  is  the 

electron  heating  rate  due  to  collisions  with  H'*',  Qf.st.o  is  the  electron  heating 
rate  due  to  fine  structure  excitation  of  neutral  oxygen,  and  Qh  is  the  heating  of 
the  electron  gas  due  to  interactions  with  photoelectrons. 

The  collisional  heating  rate  expressions  for  each  species  is  presented  below: 
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The  heating  due  to  photoelectrons  is  obtained  from  an  empirical  representa¬ 
tion  of  the  heating  rates  obtained  by  Swartz  et  al.  [1975]. 
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where  c  is  the  heating  efficiency,  A  is  the  geographic  latitude, 
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The  constzuits  are  k^  =  6,64158  x  10  k^  =  6.201753,  Jke  =  2.44626  x  10  and 
^7  =  2.1153846. 

The  final  step  is  to  convert  T,  to  a  variable  7,,  which  is  defined  as 
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and  perform  a  time-varying  Taylor’s  expansion  to  produce  a  linear  (in  terms  of 
7j),  second-order  differential  equation  of  the  form. 
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CHAPTER  IV 


RESULTS  AND  DISCUSSION 


Boundary  Conditions 


Before  the  continuity-momentum  equation  or  the  energy  equation  can  be 
solved  numerically,  appropriate  boundary  conditions  must  be  specified.  For  our 
model,  these  boundary  conditions  are  the  number  densities  for  H'*’  and  O"*",  and 
the  temperatures  of  H"^,  O'*",  and  the  electrons  at  the  ends  of  the  field  line  in  con¬ 
jugate  hemispheres.  This  is  fortuitous  since  densities  at  and  below  the  peak 
have  been  measured  and  studied  for  over  40  years  through  the  use  of  ionosondes. 
In  addition,  ion  and  electron  temperatures  in  the  F-region  have  been  measured 
using  incoherent  backscatter  radars  for  over  30  years  [Raitt,  1988].  Incoherent 
radar  facilities  at  Arecibo,  Puerto  Rico  (dip  latitude  =  30°N)  and  Jicamarca, 
Peru  (dip  latitude  =  1°N)  regularly  monitor  the  low-latitude  ionosphere  [5chunik 
and  Nagy,  1980].  At  high  latitudes,  the  field  line  can  be  assumed  to  be  nearly 
vertical,  and  the  appropriate  boimdau’y  conditions  are  the  density /temperature 
at  the  base  of  the  field  line  and  a  flux  at  the  top  of  the  field  line.  Because 
these  fluxes  are  not  well  known,  we  feel  our  model  both  includes  the  important 
physical  processes  and  utilizes  boundary  conditions  that  can  be  more  accurately 
specified. 

Our  boundary  conditions  are:  (a)  T(0''"),  T(H''‘)  =  T„  at  both  ends,  (b)  Tj 
=  1525°K,  (c)  N(0'^)  =  1  xlO^,  and  (d)  N(H'^)  =  10.  Our  values  for  T*  are 
based  on  a  synthesis  of  previous  model  results,  results  from  IRI,  and  published 
observations.  Boundary  values  are  somewhat  arbitrary  since  they  only  affect 
the  exterior  points  in  the  array. 


Inputs 

To  specify  the  neutral  atmosphere,  the  electric  field,  and  the  neutral  winds 
for  varying  levels  of  solar  activity,  several  inputs  are  required  by  the  MSIS 
model,  the  electric  field  model  [Richmond  et  al.,  1980],  and  our  neutral  wind 
model.  The  geographic  equivalents  of  the  points  along  a  given  field  line  are 
determined  using  subroutine  GGM  from  the  International  Reference  Ionosphere 
(IRI)  model.  Since  our  calculations  are  organized  in  magnetic  coordinates  and 
time  is  organized  in  geographic  coordinates,  we  assume  the  input  time  (local  or 
universal)  is  valid  for  the  northernmost  point  on  the  field  line  and  allow  time 
to  vary  along  the  field  line.  This  variation  is  in  accordance  with  the  change  in 
geographic  longitude  as  magnetic  colatitude  increases.  In  addition,  we  assume 
quiet  magnetic  conditions.  For  sunspot  minimum  conditions,  the  sunspot  num¬ 
ber  utilized  in  equation  (16)  is  zero  (Pqo  =  I-l  ^  10~^)  and  the  10.7  cm  flux  is  70; 
for  solar  maximum  conditions,  the  sunspot  number  would  be  161,  corresponding 
to  Poo  =  3.3  X  10~^,  and  the  10.7  cm  flux  is  assumed  to  be  200.  See  Figure  6 
for  results  from  MSIS86  for  sunspot  minimum  conditions  along  a  field  line  with 
an  equatorial  crossing  altitude  of  1000  km.  Figure  6a  shows  densities  for  the 
northern  hemisphere  portion  of  the  field  line  for  equinox  conditions  (Julian  day 
=  80).  The  southern  hemisphere  profile  (not  shown)  exhibits  little  variation 
from  the  results  in  the  northern  hemisphere.  Figure  66  is  valid  for  June  solstice 
(Julian  day  =  171)  conditions,  but  once  again  there  is  little  variation  between 
hemispheres  or  between  solstice  and  equinox  densities. 

Based  on  previous  calculations,  estimates  from  IRI,  and  observations,  our 
chosen  value  for  e  in  equation  (109)  is  e  =  .22.  It  should  be  emphasized  that  £  is 
multiplied  times  the  logzirithm  of  the  Qh-  Since  our  expression  for  heating  due 
to  the  photoelectrons  in  highly  empirical,  no  physical  significance  is  attached  to 
this  choice. 
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Results 


Since  the  focus  of  this  thesis  is  the  development  of  a  theoretical  model  of 
the  topside,  low-latitude  ionosphere,  presentation  of  results  will  be  confined  to 
demonstrating  the  model’s  ability  to  reproduce  physical  processes  in  a  reasonable 
manner.  It  is  assumed  that  if  the  physical  processes  are  handled  correctly, 
global  or  regional  modelling  of  the  low-latitude  ionosphere  is  a  straightforward 
application. 


Figure  7.  Dipole  field  lines  with  equatorial  crossing  altitudes  of  1000  km  and 
3000  km. 


Figure  7  shows  the  two  field  lines  being  considered.  As  previously  noted,  field 
lines  at  these  latitudes  are  inclined  to  the  vertical  and  become  horizontal  at 
the  magnetic  equator.  It  is  this  inclination  of  the  field  lines  that  contributes 
to  the  behavior  of  the  low-latitude  ionosphere.  Note  that  the  field  line  with 
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an  equatorial  altitude  of  1000  km  intersects  200  km  at  ±18°  magnetic  latitude, 
and  the  field  line  crossing  the  equator  at  3000  km  intersects  200  km  at  ±33° 
magnetic  latitude. 

As  a  test  of  the  model,  solar  minimum  conditions  for  March  equinox  and  June 
solstice  were  chosen.  Most  results  are  for  the  field  line  with  an  equatorial  crossing 
altitude  of  1000  km  (1.16  R^)  at  0°  magnetic  longitude.  However,  results  for 
a  field  line  with  an  equatorial  crossing  altitude  of  3000  km  are  also  shown  for 
comparison.  All  profiles  represent  a  steady  state  solution  of  the  model  for  local 
noon  (at  the  northernmost  point)  during  quiet  magnetic  conditions. 

Figure  8  shows  model  calculations  for  equinox  when  neutral  wind  effects, 
electrodynamic  drift,  and  adiabatic  effects  are  ignored.  The  peak  occurs 
at  an  altitude  of  about  300  km  with  a  peak  density  of  2.5  x  10®  cm“®.  O"*" 
remains  the  dominant  ion  up  to  an  altitude  of  720  km,  which  is  the  approximate 
O-H  transition  height.  The  profiles  are  fairly  symmetric  about  the  equator 
with  the  minimum  in  density  being  displaced  slightly  into  the  northern 
hemisphere  due,  perhaps,  to  the  slight  asymmetry  in  the  neutral  atmosphere. 
At  the  magnetic  equator,  densities  drop  to  5  X  10®  cm"®  and  the  density 
increases  to  1  x  10^  cm"®.  The  temperatures  of  the  ions  and  electrons  increase 
with  height  above  the  F2  peak  and  reach  maxima  of  2750°K  for  the  electrons, 
2500°K  for  the  H'*‘  ions,  and  2472°K  for  the  ions.  Below  the  F^  peak,  the  ions 
cool  to  the  neutral  temperature  due  to  the  frequent  collisions  with  the  neutrals, 
which  act  as  a  heat  sink.  The  electron  temperature  shows  a  minimum  at  the 
F2  peak  because  there  exists  a  direct  relationship  between  the  electron  density 
and  electron  temperature.  The  more  numerous  the  ambient  electron  density, 
the  less  heat  an  individual  electron  receives  from  the  photoelectrons.  This  case 
will  serve  as  a  reference  for  comparison. 


Figure  9  shows  model  calculations  for  solstice  conditions.  Although,  it’s 
not  readily  apparent,  careful  examination  reveals  that  densities  have  decreased 
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Figure  8.  Plasma  (a)  densities  and  (b)  temperatures  along  a  magnetic  field 
line  with  an  equatorial  crossing  altitude  of  1000  km  for  equinox  and  solar  mini¬ 
mum  conditions  and  0°  magnetic  longitude.  Electrodynamic,  neutral  wind,  and 
adiabatic  effects  are  neglected. 


Temperature 


i! 


w’^  'jy. 


olightly  along  the  bulk  of  the  field  line  when  compared  to  the  equinox  case.  Den¬ 
sities  in  the  southern  hemisphere  F2  peak  have  decreased  considerably  less  than 
in  the  northern  hemisphere.  H"*"  densities  remained  unchanged  in  the  north¬ 
ern  hemisphere,  but  increased  along  the  field  line  in  the  southern  hemisphere. 
Associated  with  this  decrease  in  density,  ion  and  electron  temperatures  have 
increased  in  both  hemispheres  below  600  km.  The  increase  in  ion  tempera¬ 
ture  results  directly  from  the  increase  in  electron  temperature  and  indirectly 
as  a  result  of  the  density  change.  The  transition  height  has  decreased 

about  50  km  for  solstice  conditions.  This  comparison  between  Figures  8  and 
9  demonstrate  that  the  asymmetry  between  northern  and  southern  hemisphere 
F-regions  is  not  due  to  asymmetries  in  the  neutral  atmosphere  or  differing  solar 
zenith  angles. 

Now  that  the  characteristics  of  the  reference  profiles  have  been  established,  it 
is  possible  to  evaluate  the  relative  effects  of  other  processes.  The  first  process  to 
be  considered  is  the  electrodynamic  drift,  which  is  about  2  meters/second  at  the 
equator  (They  are  assumed  to  decrease  with  distance  away  from  the  equator.). 
Results  from  our  model  with  these  effects  included  are  shown  in  Figure  10  for 
equinox  conditions  and  Figure  11  for  solstice  conditions.  Neutral  wind  effects 
and  adiabatic  effects  due  to  field-aligned  flow  are  neglected.  When  B  x  B  effects 
are  included,  the  O"*"  and  densities  are  decreased  above  340  km,  but  increase 
slightly  between  the  F2  peak  and  340  km.  Depletions  in  plasma  density  are 
greatest  at  the  magnetic  equator.  This  effect  is  consistent  with  models  of  the 
formation  of  the  Appleton  anomaly  where  plasma  drifts  upward  at  the  equator 
due  to  electrodynamic  forces  and  then  diffuses  down  the  field  lines  forming 
regions  of  enhanced  plasma  density  at  ±16°  [Anderson,  19736].  The 
transition  height  is  raised  about  65  km.  The  electron  and  ion  temperatures 
decrease  30-50°K  due  to  adiabatic  expansion  of  the  plasma. 

Figure  12  contains  profiles  of  density  and  temperature  from  the  field  line  with 
an  equatorial  crossing  altitude  of  3000  km  (1.47  R«)  at  0°  magnetic  longitude 
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Figure  9.  Plasma  (a)  densities  and  (b)  temperatures  along  a  magnetic  field 
line  with  an  equatorial  crossing  altitude  of  1000  km  for  solstice  and  solar  mini¬ 
mum  conditions  and  0°  magnetic  longitude.  Electrodynamic,  neutral  wind,  and 
adiabatic  effects  are  neglected. 
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Figure  11.  Plasma  (a)  densities  and  (b)  temperatures  along  a  magnetic  field  line 
with  an  equatorial  crossing  altitude  of  1000  km  for  solstice  and  solar  minimum 
conditions  and  0°  magnetic  longitude.  Neutral  wind  and  adiabatic  effects  are 
neglected.  Electrodynamic  drifts  included. 
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Figure  12.  Plasma  (a)  densities  and  (b)  temperatures  along  a  magnetic  field 
line  with  an  equatorial  crossing  altitude  of  3000  km  for  equinox  and  solar  mini¬ 
mum  conditions  and  0°  magnetic  longitude.  Electrodynamic,  neutral  wind,  and 
adiabatic  effects  are  neglected. 
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and  equinox  conditions.  As  with  Figure  8,  effects  due  to  electrodynamic  drift, 
neutral  winds,  and  adiabatic  effects  due  to  field-aligned  flow  are  neglected.  For 
this  field  line,  peak  F2  region  densities  are  5.2  x  10^  cm“^  and  the  height  of  the 
F2  peak  is  265  km.  At  3000  km,  the  density  drops  below  1  cm~^  and  the 
H'*"  density  is  3  x  10^  cm  The  ion  temperatures  are  nearly  equal  at  3000  km. 
The  electron  temperature  is  about  50°  warmer  than  the  ion  temperature. 

Figure  13  is  also  valid  at  solar  minimum  for  a  field  line  with  an  equatorial 
altitude  of  3000  km,  but  is  valid  for  northern  hemisphere  summer.  Once  again 
neutral  wind  effects  and  electrodynamic  drifts  are  neglected.  The  peak  den¬ 
sities  in  the  F2  region  now  exhibit  a  slight  asymmetry  between  hemispheres. 
The  peak  density  in  the  southern  hemisphere  becomes  the  larger  peak  and  the 
northern  hemisphere  peak  density  has  actually  decreased  in  magnitude  from  the 
equinox  case.  Temperatures  of  all  species  have  increased  slightly  in  the  northern 
hemisphere  and  decreased  slightly  in  the  southern  hemisphere. 

Figure  14  is  valid  for  the  same  conditions  cis  Figure  13  except  that  Figure  14 
is  valid  for  sunspot  maximum.  Most  evident  is  the  increase  in  F2  peak  densities 
from  about  5  X 10^  cm~^  to  over  1  x  10®  cm~®.  O"*"  densities  have  increased  by  two 
orders  of  magnitude  at  the  top  of  the  field  line.  H"*"  densities  increzised  slightly  at 
the  equator,  but  dropped  sharply  between  200  -  1200  km.  Temperatures  within 
the  F2  peak  drop  to  the  neutral  temperature  due  to  strong  coupling  between  the 
ions  and  the  neutrals  and  the  marked  increase  in  electron  density.  Temperatures 
of  all  species  at  the  magnetic  equator  have  increased  from  over  3000°K  to  over 
4500°K. 

Next,  the  effects  of  a  field-aligned  plasma  flow  (induced  by  a  horizontal,  merid¬ 
ional  neutral  wind)  on  the  plasma  densities  along  the  field  line  will  be  considered. 
Figure  15  shows  the  effects  of  a  summer  to  winter  hemisphere  neutral  wind  on 
plasma  densities.  As  plasma  is  pushed  up  the  field  line  in  the  northern  hemi¬ 
sphere,  it  enters  a  region  of  lower  recombination  rates  and  so  plasma  densities 
are  increaised  along  the  northern  half  of  the  field  line.  The  downward  flow  of 
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Figure  13.  Plasma  (a)  densities  and  (b)  temperacares  along  a  magnetic  field 
line  with  an  equatorial  crossing  altitude  of  3000  km  for  solstice  and  solar  mini¬ 
mum  conditions  and  0°  magnetic  longitude.  Electrodynamic,  neutral  wind,  and 
adiabatic  effects  are  neglected. 
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Figure  14.  Plasma  (a)  densities  and  (b)  temperatures  along  a  magnetic  field 
line  with  an  equatorial  crossing  altitude  of  3000  km  for  solstice  and  solar  maxi¬ 
mum  conditions  and  0°  magnetic  longitude.  Electrodynamic,  neutral  wind,  and 
ad'rbatic  effects  are  neglected. 
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Figure  15.  Plasma  densities  along  a  magnetic  field  line  with  an  equatorial  cross¬ 
ing  altitude  of  3000  km  for  solstice  and  solar  minimum  conditions  and  0°  mag¬ 
netic  longitude.  Electrodynamic  and  adiabatic  effects  are  neglected.  Neutral 
wind  effects  on  densities  are  included. 
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plasma  in  the  southern  hemisphere,  likewise,  pushes  the  plasma  into  regions  of 
higher  chemical  recombination  rates  and  plasma  densities  along  the  field  line  are 
reduced. 


CHAPTER  V 

SUMMARY  AND  CONCLUSIONS 

A  number  of  unique  transport  effects  occur  in  the  low-latitude  ionosphere  in 
the  height  range  200  -  3000  km.  These  transport  effects  include  the  formation 
of  the  Appleton  anomaly,  which  is  due  to  upward  electrodynamic  drift  of  the 
plasma  at  the  magnetic  equator  and  downward  diffusion  of  this  plasma  along 
magnetic  field  lines  during  the  day.  Because  the  field  lines  are  inclined  to  the 
vertical,  meridional  neutral  winds  can  cause  a  transport  of  ionospheric  plasma 
parallel  to  B  and  produce  asymmetries  in  the  F-region  densities  and  tempera¬ 
tures  between  conjugate  hemispheres.  Longitudinal  asymmetries  in  the  F-region 
can  arise  due  to  the  offset  in  the  magnetic  and  geographic  coordinate  systems. 
Nighttime  F-region  densities  are  maintained  at  much  higher  densities  than  can 
be  explained  by  recombination.  This  maintenance  of  the  nighttime  F-region  can 
be  explained  by  upward  transport  of  plasma  into  regions  of  lower  recombination 
rates. 

Our  model  provides  time-dependent,  self-consistent  solutions  for  the  coupled 
continuity-momentum  and  energy  equations  for  O"^,  H"*",  and  electrons  It  as¬ 
sumes  an  ionosphere  composed  of  two  primary  ion  species  (H"^  and  O’*")  and 
electrons  for  the  range  of  altitudes  under  consideration.  Additionally,  we  as¬ 
sume  the  plasma  behaves  as  a  displaced  Maxwellian  and,  therefore,  the  fluid 
equations  apply. 

The  basis  for  our  continuity-momentum  equation  is  the  work  of  Sterling  et  al. 
[1969].  Sterling  et  al.  [19691  developed  a  set  of  diffusion  coefficients  appropriate 
for  curved,  dipole  field  lines  for  a  coordinate  system  which  E  x  B  drifts  with 
the  field  line.  As  a  result,  in  this  frame  of  reference,  all  motions  are  field-aligned 
and  the  problem  becomes  one-dimensional  (parallel  to  B).  Unfortunately,  their 
model  suffers  from  overly  simplistic  assumptions  concerning  low-latitude  electric 
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fields  and  meridional,  neutral  wind  speeds.  The  energy  equation  follows  the  basic 
approach  of  Bailey  et  al.  [1973). 

The  linearized,  continuity-momentum  equation  presented  couples  calcula¬ 
tions  of  H"^  and  densities  and  diffusion  velocities.  To  solve  the  continuity- 
momentum  equation  and  the  energy  equations,  we  adopt  an  orthogonal,  curvi- 
linezir  coordinate  system  in  which  the  dipole  field  lines  are  lines  of  constant  p  and 
the  potentials  are  lines  of  constant  q.  Any  point  can  now  be  uniquely  defined 
in  terms  of  p  and  q.  If  we  calculate  along  a  field  line  (p  =  constant),  then  all 
variations  can  be  expressed  in  terms  of  q. 

Dominant  chemistry  in  the  200  -  3000  km  height  range  is  included.  These 
reactions  involve  the  photoionization  of  neutral  O  and  the  charge  exchange  re¬ 
actions  between  and  Although  only  important  at  the  lowest  heights 
being  considered,  the  loss  of  in  reactions  with  N2  and  O2  are  also  included. 
Calculations  of  the  photoionization  of  O  use  the  Chapman  grazing  incidence 
function  [Smith  and  Smith,  1972). 

Flow  velocities  zure  obtained  from  the  momentum  equation.  To  solve  for  these 
flow  velocities,  we  assume  steady  state,  that  viscous  and  nonlinear  terms  can  be 
neglected,  and  that  all  flow  is  parallel  to  B.  The  electric  field  in  the  momentum 
equation  for  each  ion  species  is  calculated  from  the  electron  momentum  equa¬ 
tion  by  assuming  the  flow  of  the  electrons  is  determined  by  a  balance  between 
electrostatic  and  pressure  gradient  forces.  The  appropriate  collision  frequencies 
are  obtained  from  Schunk  [1983).  These  flow  velocities  are  then  substituted  into 
the  continuity  equation. 

Next,  the  continuity  equation  is  evaluated  to  yield  ion  densities  along  the 
field  line.  The  velocity  of  each  ion  species  is  separated  into  its  parallel  and 
perpendicular  components.  The  perpendicular  component  of  drift  is  incorpo¬ 
rated  with  the  local  time  derivative  to  yield  the  convective  derivative  (fixing  the 
coordinate  system  on  the  field  line).  Divergences  and  diffusion  coefficients  are 
then  evaluated  using  Sterling  et  al's  [1969)  methods,  except  we  chose  to  finite 
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difference  some  derivatives  as  opposed  to  expanding  each  term  (as  Sterling  et 
al.  [1969]  did).  Once  these  divergences  are  substituted  back  into  the  continuity 
equation,  a  nonlinear,  second-order  partial  differential  equation  results.  To  solve 
this  equation  numerically,  it  must  first  be  linearized.  This  is  accomplished  by 
performing  a  time-varying  Taylor’s  expansion  on  all  nonlinear  terms. 

The  methods  and  assumptions  used  in  deriving  a  linearized  energy  equation 
for  each  species  are  much  the  same  as  those  for  the  continuity-momentum  equa¬ 
tion.  Our  derivation  begins  with  the  general  form  of  the  energy  equation  found  in 
Sehunk  [1983].  Using  the  assumption  that  the  plasma  approximates  a  displaced 
Maxwellian,  the  5-moment  approach  can  be  used  to  calculate  the  rate  of  energy 
exchange  due  to  collisions.  Next,  the  form  of  the  continuity  equation  valid  in  a 
diffusive  region  is  used  in  conjunction  with  the  definition  of  the  heat  flow  vector 
to  obtain  the  time  rate  of  change  of  the  temperature  of  a  given  species.  Within 
the  energy  equation,  account  is  taken  of  the  frictional  heating/cooling  of  the  ion 
species  due  to  relative  motion  between  the  ion  gas  and  the  neutral  gas,  adiabatic 
expansion  and  compression,  conduction  along  the  field  lines,  and  collisional  heat 
exchange.  In  general,  the  O’*"  ions  gain  heat  through  collisions  with  the  electron 
and  H"*"  gases,  and  lose  heat  through  inelastic  interactions  (charge  exchange  re¬ 
actions)  with  H  and  elcistic  collisions  with  the  neutrals.  The  H"*"  ions  gain  heat 
from  collisions  with  the  electrons,  and  lose  heat  through  collisions  with  O"^,  the 
neutrals,  and  O.  The  ambient  electron  gas  loses  heat  through  collisions  with 
O"^,  H‘’’,and  the  fine  structure  excitation  of  O. 

Bailey  et  al.  [1973]  only  consider  the  nighttime  case.  Since  we  were  interested 
in  modelling  a  variety  of  local  times,  a  heating  rate  due  to  interaction  of  the 
photoelectrons  with  the  ambient  electron  gas  had  to  be  adopted  for  daytime 
conditions.  An  empirical  relation  was  developed  from  the  results  of  Swartz  et 
al.  [1975].  Initially,  the  steady  state  electron  temperatures  obtained  using  this 
empirical  model  were  much  higher  than  previous  calculations  and  observations. 
The  heating  factors  were,  therefore,  scaled  to  obtain  results  more  consistent 
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previously  reported  values. 

Once  the  divergences  are  evaluated  in  the  energy  equations,  a  second-order, 
partial  differential  equation  for  each  species  results.  Linearization  is  accom¬ 
plished  by  converting  the  temperatures  to  an  independent  variable,  7.  The 
equations  are  then  linearized  in  terms  of  7  using  time-varying  Taylor’s  expan¬ 
sion. 

A  chief  criticism  of  previous  modelling  efforts  is  that  they  employ  very  sim¬ 
plistic  electric  fields  or  neutral  wind  distributions.  While  these  techniques  may 
produce  more  accurate  results  in  an  isolated  case,  they  are  useless  in  general 
modelling  applications  where  large  errors  would  certainly  occur.  We  attempt 
to  improve  upon  previous  efforts  by  incorporating  average  neutral  wind  speeds 
obtained  from  the  Thermospheric  Global  Circulation  Model  (TGCM)  [Roble  et 
al.,  1977]  for  solar  maximum  and  minimum  conditions  at  equinox  and  solstice. 
Daily  variations  were  introduced  by  interpolating  wind  speeds  between  equinox 
and  solstice  conditions  based  on  the  solar  declination  angle.  When  mapped 
into  geomagnetic  coordinates,  longitudinal  variations  in  the  mean  neutral  wind 
speeds  axe  introduced  because  of  the  tilt  in  the  dipole  field  relative  to  the  earth’s 
rotational  axis. 

To  estimate  electric  fields  responsible  for  the  electrodynamic  drift  of  iono¬ 
spheric  plasma,  we  adopt  the  empirical  model  of  Richmond  et  al  [1980].  Devel¬ 
oped  for  solar  minimum  conditions,  the  incorporation  of  this  model  allows  for 
much  more  realistic  estimates  of  the  electric  fields  and  the  electrodynamic  drifts 
than  previously  possible.  With  the  help  of  Richmond  [1988],  a  method  was  de¬ 
vised  to  calculate  variations  of  this  electrodynamic  drift  with  height.  Previous 
authors  [Anderson,  1973a,  19736;  Chan  and  Walker,  1984a,  19846]  had  assumed 
an  electric  field  or  electrodynamic  drift  velocity,  and  assumed  it  was  independent 
of  altitude.  Our  calculations  for  solar  minimum,  June  solstice,  0°  magnetic  lon¬ 
gitude,  and  local  noon  show  that  the  variation  of  the  drift  velocity  with  height 
for  the  field  line  crossing  the  magnetic  equator  at  1000  km  is  approximately 


one-half  the  vertical  drift  velocity.  Since  they  are  of  the  same  order  of  magni¬ 
tude,  it  is  not  clear  that  neglecting  this  variation  of  vertical  drift  velocity  with 
altitude  is  a  valid  assumption.  The  east-west  drift  has  generally  been  ignored 
based  on  the  results  of  Sterling  et  al.  (1969];  but  since  his  electric  fields  were 
very  simplistic,  the  importance  of  east-west  electrodynamic  drifts  needs  to  be 
re-evaluated.  Since  we  ran  only  cases  for  local  noon,  evaluation  of  the  east-west 
electrodynamic  drift  will  be  undertaken  in  future  tests  of  the  model. 

The  prime  focus  of  this  research  is  the  development  of  a  low-latitude  iono¬ 
spheric  model.  Results  were,  therefore,  aimed  at  demonstrating  that  the  model 
could  reproduce  physical  processes  in  the  low-latitude  ionosphere.  Steady  state 
profiles  of  density  and  temperature  were  created  for  the  magnetic  prime  merid¬ 
ian  for  solar  minimum,  March  equinox,  and  local  noon.  Electrodynamic,  neutral 
wind,  and  adiabatic  effects  were  not  included  in  the  model  computations.  This 
case  served  as  a  reference  for  the  test  cases  to  follow.  When  this  case  was 
repeated  for  solstice  conditions,  very  little  asymmetry  was  exhibited  between 
hemispheres,  but  then  there  was  little  asymmetry  in  the  neutral  atmosphere. 
This  test  shows  that  the  asymmetry  between  the  F-region  densities  that  has 
been  observed  during  solstice  conditions  is  not  solely  due  to  asymmetric  solar 
input. 

A  second  test  involved  the  inclusion  of  electrodynamic  drifts  for  equinox  and 
solstice  conditions.  Density  profiles  showed  the  expected  depletions  of  pleisma  at 
the  magnetic  equator  and  an  increase  in  plasma  density  at  ±  18°  near  300  km. 
Plasma  temperatures  decreased  along  the  field  line  in  response  to  the  adiabatic 
expansion  of  the  plasma  as  it  was  being  forced  upwards.  Therefore,  expected 
adiabatic  effects  were  also  reproduced  correctly. 

When  strong  north-to-south  neutral  winds  were  introduced,  the  model  once 
again  reproduced  the  correct  trends.  This  effect  is  consistent  with  increased 
transport  of  plasma  to  greater  heights  where  recombination  rates  are  lowei .  For 
the  southern  hemisphere,  the  winds  resulted  in  lower  plasma  densities  because 
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more  plasma  was  being  forced  into  regions  of  higher  recombination  rates. 

In  summary,  our  model  seems  to  reproduce  the  important  physical  processes  in 
the  low-latitude  ionosphere  in  a  manner  consistent  with  theory  and  observations. 
We  have  achieved  significant  improvements  over  previous  modelling  efforts  by 
including  more  realistic  estimates  of  neutral  winds  and  electric  fields.  In  the 
process,  we  have  increased  the  general  applicability  of  our  model  to  a  wide 
variety  of  seasons,  levels  of  solar  activity,  longitudes,  and  local  times.  An  area 
of  further  improvement  for  this  model  would  be  to  include  a  more  physical 
treatment  of  the  photoelectrons. 
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APPENDICES 


Appendix  A:  Neutral  Wind 
Model 


The  neutral  wind  model  used  in  this  study  allows  for  variations  of  the  north- 
south  component  of  the  neutral  wind  speed  with  magnetic  latitude,  magnetic 
longitude,  altitude,  day  of  the  year,  and  local  time.  To  accomplish  this,  the 
analytic  expressions  for  the  diurnal,  V24  (Al),  semi-diurnal,  V12  (A2),  terdi- 
urnal,  V8  (A3)  variations,  and  the  total  wind  speed,  U  (A5),  were  taken  from 
Anderson  et  al.,  [1987]  except  that  $24  (A4)  has  been  redefined. 


V24  =  2^24  cos  -  18  - 

(Al) 

K12  =  50.m('  3  '“')1 

V180  3  y  [12  V  12  J  _ 

(A2) 

r  27r 

V8  =  0.38a cos  I (t  -  8.2  +  0.079  |a|) 

(A3) 

1^1 

=  “+'0,^3  5 

(Al) 

i7  =  K0-f  K8  +  712  +  T24 

(AS) 

S  =  arcsin[23.5sin(i2)j 

(A6) 

(A7) 

"“(w) 

(AS) 

where  a  is  the  geographic  latitude,  t  is  the  local  time,  6  is  the  solar  declination 
angle,  DD  is  the  Julian  date,  R  is  the  number  of  days  after  March  equinox 
converted  to  degrees,  and  U  is  the  resulting  meridional  wind  speed  in  meters 
per  second  defined  positive  from  north  to  south. 
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The  mean  wind  speed,  VO,  was  obtained  from  the  results  of  Roble  et  al.  [1977] 
in  which  they  give  mean  values  for  the  zonal  and  meridional  components  of  the 
neutral  wind  for  equinox  and  solstice  conditions  for  solar  minimum  and  solar 
maximum  in  geographic  coordinates.  Estimates  of  these  values  for  every  10° 
latitude  from  40°  5  to  AQ°N  and  every  50  km  from  200  km  to  600  km  are  shown 
in  Tables  1-6.  Since  the  results  of  Roble  et  al.  [1977]  did  not  extend  to  600 
km  in  all  cases,  winds  are  assumed  to  be  constant  above  the  maximum  height 
presented  in  their  paper.  This  maximum  height  varies  from  400  km  for  solstice, 
solar  minimum  conditions  to  600  km  for  equinox,  solar  maximum  conditions. 
For  points  outside  this  range  of  altitudes  or  geographic  latitudes,  the  winds  are 
assumed  constant  and  equal  to  the  value  for  the  mean  wind  at  the  nearest  point 
inside  the  array.  Since  our  model  requires  the  component  of  the  neutral  wind 
speed  along  the  magnetic  field  line,  the  components  of  the  neutral  wind  speed  in 
geographic  coordinates  are  converted  to  magnetic  coordinates  using  equations 
(A9)  and  (AlO)  given  below. 

Ueq  =  U^onaU  shl  (A)  +  Umeu  COS  (A)  (A9) 

Utol  -  U zonal,  sin  (A)  +  Umer,  COS  (A)  (A  10) 

where  A  is  the  magnetic  declination  angle,  Ugonal,  is  the  zonal  component  of  the 
neutral  wind  speed  at  equinox,  Umer,  is  the  meridional  component  of  neutral 
wind  speed  at  equinox,  U zonal,  is  the  zonal  component  of  the  neutral  wind  speed 
at  solstice,  and  Umer,  is  the  meridional  component  of  the  neutral  wind  speed 
at  solstice.  Now,  the  mean,  meridional  component  of  the  neutral  wind  speed, 
VO  (All),  can  be  determined  by  interpolating  between  the  solstice  and  equinox 
values  based  on  the  solar  declination  angle. 

VO  =  Uz,  ~  (^(Uz,  -  U,oi)  (^11) 
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TABLE  1.  Zonal  Neutral  Wind  Speeds  for  Equinox 
and  Solar  Maximum  Conditions.  Speeds  are  in 
Meters/Second  and  Eastward  Winds  are 
Defined  as  Positive. 
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TABLE  2.  Meridional  Neutral  Wind  Speeds  for  Equinox 
and  Solar  Maximum  Conditions.  Speeds  are  in 
Meters/Second  and  Southward  Winds  are 
Defined  as  Positive. 
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TABLE  3.  Zonal  Neutral  Wind  Speeds  for  Winter 
Solstice  and  Solar  Maximum  Conditions. 
Speeds  are  in  Meters/ Second  and 
Eastward  Winds  are  Defined  as 
Positive 
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TABLE  4.  Meridional  Neutral  Wind  Speeds  for  Winter 
Solstice  and  Solar  Maximum  Conditions.  Speeds  are 
in  Meters/Second  and  Southward  Winds  are  Defined 
as  Positive. 
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Finally,  the  total  meridional  component  of  the  neutral  wind  speed  is  the  sum 
of  the  mean  wind  speed  (All),  the  diurnal  variation  (Al),  the  semi-diurnal 
variation  (A2),  and  the  terdiurnal  variation  (A3)  as  given  by  equation  (AS). 

The  results  for  a  variety  of  input  conditions  are  shown  in  Figures  16  -  22. 
Figure  16  shows  the  variation  of  the  north-south  component  of  the  neutral  wind 
for  solar  maximum,  equinox  conditions  along  97°  geomagnetic  longitude  at  mid¬ 
night.  Since  magnetic  latitudes  and  geographic  latitudes  approximately  coincide 
at  this  magnetic  longitude,  the  winds  along  this  magnetic  meridian  are  fairly 
symmetric  about  the  equator  as  might  be  expected  for  equinox.  Winds  in  both 
magnetic  hemispheres  are  directed  toward  the  magnetic  equator  with  maximum 
winds  exceeding  90  meters /second  at  ±30°  magnetic  latitude,  and  are  nearly 
constant  with  height. 

Figure  17  shows  winds  for  local  noon,  but  is  otherwise  valid  for  the  same  con¬ 
ditions  as  Figure  16.  As  in  Figure  16,  the  wind  pattern  is  fairly  symmetric  about 
the  equator.  However,  unlike  Figure  16,  there  is  now  considerable  variation  of 
wind  speed  with  height  up  to  500  km  (Winds  are  assumed  constaint  above  600 
km.),  and  maximum  wind  speeds  in  the  northern  hemisphere  are  approximately 
double  those  in  the  southern  hemisphere. 

Figure  18  is  an  example  of  how  neutral  winds  can  vary  with  geomagnetic  lon¬ 
gitude  due  to  the  offset  in  the  geographic  and  geomagnetic  poles.  Conditions 
in  this  case  are  the  same  as  in  Figure  16  except  the  magnetic  longitude  is  now 
179.8° J?.  This  corresponds  to  the  magnetic  longitude  where  the  magnetic  lati¬ 
tudes  and  geographic  latitudes  are  offset  by  about  11.4°.  This  offset  results  in 
a  shifting  of  the  entire  wind  pattern  11.4°  towards  the  south  and,  hence,  the 
symmetry  between  hemispheres  is  lost.  For  this  case,  winds  do  not  reverse  until 
11.4°S  and  maximum  wind  speeds  in  the  northern  hemisphere  reach  100  me¬ 
ters/second;  whereas,  maximum  wind  speeds  in  the  southern  hemisphere  only 
reach  50  meters/second.  Therefore,  the  variation  in  wind  speed  with  magnetic 
longitude  for  a  fixed  local  time  in  solely  determined  by  the  offset  in  the  two 
coordinate  systems. 

When  compared  with  Figure  16,  Figure  19  gives  an  example  of  seasonal  vari¬ 
ations  of  neutral  wind  speeds  along  a  magnetic  field  line.  Figure  19  is  valid 
for  essentially  the  same  conditions  as  Figure  16  except  the  season  is  now  June 
solstice.  For  this  case,  winds  are  generally  directed  from  the  summer  hemi¬ 
sphere  to  the  winter  hemisphere,  but  they  do  reverse  near  15°S.  Wind  speeds 


Figure  16.  Meridional  component  of  the  neutral  winds  in  meters/second  along 
a  magnetic  field  line.  Winds  are  positive  from  north  to  south  and  the  contour 
interval  is  20  meters/second.  The  winds  are  valid  for  97° E  magnetic  longitude, 
equinox,  solar  maximum,  and  midnight  local  time. 
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Figure  17.  Meridional  component  of  the  neutral  winds  in  meters /second  along 
a  magnetic  field  line.  Same  as  Figure  16  except  the  local  time  is  1200.  Winds 
are  positive  from  north  to  south. 
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Figure  18.  Meridional  component  of  the  neutral  winds  in  meters/second  along  a 
magnetic  field  line.  Winds  are  positive  from  north  to  south  and  the  contour  in¬ 
terval  is  20  meters/second.  The  winds  are  valid  for  179.8°E  magnetic  longitude, 
equinox,  solar  maximum,  and  midnight  local  time. 
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Figure  19.  Meridional  component  of  the  neutral  winds  in  meters/second  along 
a  magnetic  field  line.  Winds  are  positive  from  north  to  south.  The  winds  are 
valid  for  91° E  magnetic  longitude,  solar  maximum,  June  solstice,  and  midnight 
local  time.  Contour  intervals  are  20  meters/second. 


have  increased  considerably  from  the  Equinox  case  in  the  northern  hemisphere 
with  a  mzocimum  in  excess  of  170  meters/second,  and  decreased  to  about  60 
meters/second  in  the  southern  hemisphere.  For  the  midnight  case  at  solstice, 
the  winds  increase  with  altitude  up  to  500  km  above  which  they  are  assumed  to 
be  constant,  and  show  considerable  asymmetry  as  would  be  expected. 

Figure  20  gives  some  typical  neutral  winds  at  1200  LT  for  the  same  conditions 
as  Figure  19.  Winds  in  this  case  are  directed  from  the  summer  to  the  winter 
hemisphere  with  the  maximum  in  the  wind  speeds  above  300  km  occurring  at  the 
latitude  corresponding  to  the  subsolar  point.  Therefore,  wind  speeds  generally 
decrease  away  from  the  latitude  of  the  subsolar  point  at  these  altitudes.  Below 
300  km,  the  winds  are  nearly  constant  with  latitude.  The  increase  in  wind  speed 
with  altitiude  at  a  given  latitude  is  maintained  from  the  midnight  case  (Winds 
are  assumed  constant  above  500  km).  Wind  speeds  for  the  local  noon  case  at 
solstice  are  much  stronger  than  for  the  equinox  case  shown  in  Figure  16,  and 
the  region  of  maximum  winds  has  shifted  from  the  lower  levels  for  equinox  to 
the  upper  levels  for  solstice. 

Figure  21  shows  estimated  wind  speeds  for  June  solstice,  solar  maximum, 
179.8° JE,  and  midnight  local  time.  This  figure  shows  the  effects  of  the  offset 
of  the  magnetic  and  geographic  axes  for  the  solstice  case.  Once  again,  the 
entire  wind  pattern  has  shifted  in  accordance  with  the  offset  between  the  two 
coordinate  systems  for  the  chosen  magnetic  longitude.  M20cimum  wind  speeds 
for  the  northern  hemisphere  have  increased  to  over  200  meters/second  amd  the 
winds  are  now  almost  exclusively  from  the  summer  to  winter  hemisphere.  The 
increase  in  wind  speeds  with  alitiude  have  become  more  pronounced  for  the 
summer  hemisphere. 

Figure  22  gives  estimated  winds  for  solar  minimum,  June  solstice,  358.6° £? 
magnetic  longitude,  and  midnight  local  time.  At  358.6° £?,  the  overall  wind 
pattern  is  shifted  about  11.4°  northward  due  to  the  offset  in  the  coordinate 
systems,  and  is,  therefore,  opposite  to  the  case  presented  in  Figure  21  with  the 
added  effect  of  decreased  solar  activity.  One  of  the  most  noticeable  results  is 
the  nearly  synunetric  pattern  the  winds  now  exhibit  which  would  normally  be 
associated  with  equinox  type  conditions  if  the  magnetic  and  geographic  poles 
were  assumed  to  coincide.  When  the  effects  due  to  shifting  the  wind  pattern 
are  removed  (by  comparing  with  Figure  19),  the  effect  due  to  decreasing  solar 
activity  on  the  wind  speeds  is  to  decrease  them  about  -  30  meters/second. 
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Figure  22.  Meridional  component  of  the  neutral  winds  in  meters/second  along 
a  magnetic  field  line.  Winds  are  positive  from  north  to  south.  The  winds 
are  valid  for  358.6°/^  magnetic  longitude,  solar  maximum,  June  solstice,  and 
midnight  local  time.  Contour  intervals  are  20  meters/second. 

Since  the  induced,  parallel  (to  B  )  drift  of  the  ions  due  to  the  neutral  winds 
can  have  significant  effects  on  ionospheric  densities  by  moving  ions  into  regions  of 
faster  or  slower  recombination  or  on  ion  temperatures  due  to  adiabatic  heating  or 
cooling,  accurately  modelling  seasonal,  local  time,  soIm  cycle,  and  longitudinal 
variations  of  the  neutral  winds  becomes  critical.  Although  this  model  only  gives 
crude  estimates  of  these  variations,  it  will  extend  previous  modelling  efforts 
[Bailey  et  al.  1973,  Bailey  and  Heelis,  1980]  over  a  wider  range  of  conditions. 


Appendix  B:  An  Empirical  Model  of 
the  Magnetic  Field 


The  magnetic  field  model,  FIELDG,  incorporated  into  this  study  was  taken 
from  the  IRI-86  empirical  model  of  the  ionosphere.  According  to  the  documen¬ 
tation  accompanying  the  subroutine  in  IRI,  FIELDG  is  a  special  version  of  the 
POGO  68/10  magnetic  field,  Legendre  model.  FIELDG  uses  transformation 
coefficients  valid  for  1973  {Sheikh  et  al.,  1978). 

Figure  23  shows  the  magnetic  field  strength  at  200  km  using  FIELDG. 
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Figure  23.  Global  magnetic  field  strengths  from  FIELDG  for  200  km.  The 
contour  labels  are  10~^  Teslas  and  the  contour  interval  is  5  x  10”®  Teslas. 

FIELDG  reproduces  the  South  Atlantic  anomaly  and  the  relative  maximum 
in  field  strength  near  70° iV,  110°i?.  Estimating  from  Figure  23,  the  magnetic 
poles  appear  to  be  located  near  60° iV,  210° E  and  60°S,  140°E. 


Appendix  C:  An  Empirical  Model  of 
Electric  Fields  at  Low  Latitudes 


Determination  of  the  electrodynamic  drift  velocities  along  magnetic  field  lines 
requires  the  electric  field  be  specified.  Therefore,  the  empirical  electric  field 
model  of  Richmond  et  al.  (1980)  was  adopted.  This  model  gives  values  for 
electrodynamic  drift  and  the  “pseudo  electrostatic  potential”  (as  defined  by 
Richmond  et  al,  [1980])  at  300  km  for  ±65°  magnetic  latitude,  universal  time, 
and  day  of  the  year.  The  model  was  constructed  using  data  from  incoherent 
radaxs  located  at  Millstone  Hill  (42.6° JV,  288.5° J?)  ,  St.  Santin  (44. 1°^,  2°E), 
Arecibo  (18.5° TV,  293.2° fT),  and  Jicamarca  (12° 5,  283.1°£T)  for  quiet  magnetic 
conditions  during  solar  minimum.  These  pseudopotentials  are  assumed  to  be  es¬ 
sentially  constant  along  the  field  lines.  This  final  point  is  important  to  our  study 
because  it  allows  estimation  of  the  electrodynamic  drift  velocities  for  heights 
other  than  300  km. 


The  estimation  process  begins  by  assuming  the  electro dynaunic  drifts  are  con¬ 
stant  between  200  km  and  300  km  because  for  field  lines  crossing  the  magnetic 
equator  below  300  km,  pseudopotentials  cannot  be  extrapolated  from  model 
calculations.  For  points  above  300  km,  the  electrodynamic  drift  velocities  are 
calculated  by  assuming  the  magnetic  field  lines  are  equipotential  and  by  using 
the  equations  for  electrodynamic  drift.  These  equations  are 
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where  $  is  the  pseudo  electrostatic  potential, 

Vx  is  the  vertical  component  of  the  electrodynamic  drift, 
is  the  azimuthal  component  of  the  electrodynamic  drift, 
I  is  the  inclination  angle, 

0  is  the  colatitude,  and 
B  is  the  magnetic  ficM  strength. 
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Equation  (Cl)  is  taken  directly  from  Richmond  et  al.,  (I980j;  however,  equa¬ 
tion  (C2)  had  to  be  modified  to  take  into  account  the  total  electric  field  since 
^  goes  to  zero  at  the  magnetic  equator  [Richmond,  1988].  Using  these  assump¬ 
tions,  the  variation  of  the  pseudopotential  with  magnetic  longitude,  altitude, 
and  colatitude  can  be  determined;  and  since  the  field  lines  are  assumed  to  be 
of  constant  potential,  the  variation  of  pseudopotential  with  longitude  is  also 
constant.  Therefore,  to  calculate  the  vertical  drift  it  only  becomes  necessary  to 
find  the  magnetic  latitude  where  the  field  line  intersects  300  km  to  determine 
along  the  entire  field  line.  Once  this  value  is  obtained,  the  vertical  drift  at 
any  point  along  the  field  line  can  be  calculated  from  equation  (Cl). 

To  obtain  values  for  the  eastward  electrodynamic  drift,  the  variation  of  pseu¬ 
dopotential  with  altitude  and  colatitude  must  be  determined.  Figure  24  shows 
a  schematic  of  how  the  variation  with  colatitude  can  be  estimated  for  any  point, 
P,  along  the  field  line.  The  process  begins  by  holding  the  radius  constant  and 
locating  the  field  lines  that  lie  from  point  P.  These  points  of  intersection 
are  indicated  as  points  A  and  B  in  Figure  24.  Once  determined,  the  field  lines 
are  traced  down  to  their  points,  C  and  D,  of  intersection  with  300  km.  Using 
the  electric  field  model,  the  pseudopotentials  at  points  C  and  D  can  next  be 
determined  and  the  variation  in  pseudopotential  can  be  obtained  by  dividing 
the  difference  in  pseudopotential  by  the  chosen  value  for  AO.  The  variation 
of  pseudopotential  is  determined  in  a  like  manner.  The  magnitude  of  the  2lz- 
imuthal  drift  is  now  computed  from  equation  (C2).  Since  the  magnitude  of  the 
electrodynamic  drifts  is  a  function  of  the  magnetic  field  strength  at  any  given 
point,  B  is  calculated  as  a  function  of  latitude,  longitude,  and  altitude  by  the 
subroutine  FIELDG  described  in  Appendix  B.  Finally,  although  the  model  al¬ 
lows  for  various  averages  in  universal  time  and  season,  none  were  used  in  our 
calculations. 

Figure  25  shows  the  vertical  electrodynamic  drift  at  300  km  as  a  function  of 
magnetic  latitude  and  universal  time  for  equinox  and  0°  magnetic  longitude. 
For  these  conditions,  the  vertical  drift  is  downward  for  virtually  all  magnetic 
latitudes  from  0000  UT  until  1000  UT,  and  is  upward  from  1000  UT  to  1800 
UT.  After  1800  UT,  the  vertical  drifts  are  upward  within  10°  of  the  equator  and 
downward  elsewhere.  Upward  drifts  peak  near  1500  UT  with  a  velocity  of  20 
meters/second;  the  downward  velocities  peak  near  0600  UT  with  a  velocity  of 
16  meters/second. 
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MAGNETIC  COLATITUDE 

Figure  24.  Schematic  representation  of  how  the  variation  of  pseudopotential  with 
colatitude  is  determined  for  any  point,  P,  along  the  magnetic  field  line.  A0  is 
defined  at  point,  P  and  A$  is  calculated  as  the  difference  in  the  pseudopotentials 
given  by  the  model  at  points,  C  and  D. 

Figure  26  is  valid  for  the  same  conditions  as  Figure  25  except  that  the  drifts 
are  now  for  1000  km.  The  pattern  of  drifts  is  essentially  the  same  as  for  the 
300  km  case  except  that  the  magnitude  has  been  reduced  about  25%.  This  is 
to  be  expected  owing  to  the  method  in  which  the  variations  with  altitude  are 
calculated.  The  method  essentially  brings  the  winds  valid  at  300  km  for  a  given 
field  line  up  to  the  altitude  of  the  point  in  question  and  reduces  them  in  an 
amount  given  by  the  variations  of  B  and  r. 

Figure  27  is  valid  for  June  solstice  at  300  km;  all  other  conditions  remain  the 
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Figure  25.  Vertical  component  of  the  electrodynamic  winds  in  meters/second  at 
300  km  for  equinox  and  0°  magnetic  longitude  as  a  function  of  magnetic  latitude 
and  universal  time.  Contours  intervals  are  20  meters/second.  Upward  winds  are 
defined  to  be  positive. 


Universal  Time 


Figure  26.  Vertical  component  of  the  electrodynamic  winds  in  meters/second 
at  1000  km  for  equinox  and  0°  magnetic  longitude  as  a  function  of  magnetic 
latitude  and  universal  time.  Contours  intervals  are  20  meters/second.  Upward 
winds  are  defined  to  be  positive. 
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same.  Peak  drift  velocities  are  somewhat  reduced  from  the  equinox  case  and  an 
eastward  bulge  is  introduced  into  the  overall  drift  pattern.  Near  the  magnetic 
equator,  the  drift  is  downward  from  2300  UT  until  1200  UT;  but  near  ±30° 
magnetic  latitude  the  drifts  shift  directions  several  times  during  the  course  of  a 
day. 


Figure  28  shows  the  vertical  drift  velocities  for  June  solstice  conditions  at  1000 
km.  As  in  Figure  26,  the  pattern  of  drifts  changes  little  with  altitudes  and  the 
magnitudes  of  the  drifts  are  reduced  about  25%. 


Even  though  of  greater  magnitude,  azimuthal  drifts  are  generally  of  much  less 
importance  and  have  been  ignored  in  calculations  of  diffusion  along  magnetic 
field  lines  [Sterling  et  al.,  1969].  However,  near  sunrise  and  sunset  horizontal 
variations  in  the  ionosphere  do  become  significant  [jRichmond,  1988]  amd  the 
azimuthal  component  is,  therefore,  retained  in  this  study. 


Figure  29  shows  the  east-west  component  of  the  electrodynaunic  drift  at  300 
km  for  equinox  and  0°  magnetic  longitude.  Westward  winds  peak  near  1600  UT 
at  50  meters/second.  There  is  an  eaistward  bulge  in  the  wind  pattern  that  is 
centered  on  the  magnetic  equator;  perhaps,  due  to  solar  zenith  angle  effects. 


Figure  30  is  valid  for  the  same  conditions  as  Figure  29  except  that  the  altitude 
is  now  1000  km.  As  in  the  case  of  the  vertical  winds,  the  east-west  winds 
essentially  exhibit  the  same  pattern  at  300  km  and  1000  km  and  the  magnitude 
of  the  winds  decrease  20  -  30%  from  300  km  to  1000  km.  Peak  winds  at  1000 
km  are  60  meters/second  in  the  eastward  direction  and  30  meters /second  in  the 
westward  direction.  Variations  of  the  wind  speed  and  direction  with  universal 
time  are  nearly  the  same  as  those  at  300  km  except  as  already  noted. 


Figure  31  shows  the  azimuthal  component  of  the  electrodynamic  drift  veloci¬ 
ties  at  300  km  for  June  solstice  2md  0°  magnetic  longitude.  For  the  solstice  case, 
winds  are  predominantly  eastward  along  this  magnetic  meridian  and  reach  a 
peak  of  over  80  meters/second  at  0200  UT.  The  drift  is  westward  only  between 
1000  and  1800  UT.  Peak  westward  winds  are  40  meters/second  and  occur  near 
1500  UT. 


Figure  32  is  the  same  as  Figure  31  except  the  winds  are  valid  for  1000  km 
versus  300  km.  Once  again,  windspeeds  decrease  20  -  30%  from  300  to  1000 
km;  and  because  of  this  decrease,  westward  drifts  in  excess  of  20  meters/second 
are  confined  to  a  band  of  only  3  houre  in  universal  time.  At  1000  km,  the  peak 
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Figure  27.  Vertical  component  of  the  electrodynaunic  winds  in  meters/second 
at  300  km  for  June  solstice  and  0“  magnetic  longitude  as  a  function  of  magnetic 
latitude  and  universal  time.  Contours  intervals  are  20  meters/second.  Upward 
winds  are  defined  to  be  positive. 
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Figure  28.  Vertical  component  of  the  electrodynamic  winds  in  meters/second 
at  1000  km  for  June  solstice  and  0°  magnetic  longitude  as  a  function  of  magnetic 
latitude  and  universal  time.  Contours  intervals  are  20  meters/second.  Upward 
winds  are  defined  to  be  positive. 
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Figure  29.  Azimuthal  component  of  the  electrodynamic  winds  in  meters/second 
at  300  km  for  equinox  and  0®  magnetic  longitude  as  a  function  of  magnetic 
latitude  and  universal  time.  Contours  intervals  are  20  meters/second.  Eastward 
winds  are  defined  to  be  positive. 
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Figure  30.  Azimuthal  component  of  the  electrodynamic  winds  in 
meters/second  at  1000  km  for  equinox  and  0®  magnetic  longitude  as  a  function  of 
magnetic  latitude  and  universal  time.  Contours  intervals  are  20  meters/second. 
Eastward  winds  are  defined  to  be  positive. 


westward  winds  only  reach  about  2f>  meters/second,  whereas,  the  eastward  drift 
exceeds  65  meters/second. 

In  summary,  azimuthal  wind  speeds  are  considerably  greater  than  the  ver¬ 
tical  drifts,  but  owing  to  the  nearly  horizontally  homogeneous  nature  of  the 
ionosphere  throughout  much  of  the  day,  they  are  generally  less  signihcant  to 
the  consideration  of  diffusion  of  ionospheric  plasma  along  magnetic  field  lines. 
Upward  vertical  drifts  tend  to  be  greater  for  equinox  than  for  June  solstice, 
whereas,  downward  drifts  vary  little  with  season  for  the  case  considered.  In  the 
case  of  the  azimuthal  drifts,  this  pattern  is  reversed  with  the  faster  drifts  oc¬ 
curring  for  solstice  rather  than  equinox  conditions.  While  winds  are  symmetric 
about  the  magnetic  equator  in  all  cases,  the  vertical  drift  component  exhibits 
more  latitudinal  variation  than  the  azimuthal  component.  Finally,  it  is  worth 
noting  again  that  the  cases  presented  are  meant  to  be  representative  of  quiet 
magnetic  conditions  and  low  solar  activity  along  the  magnetic  prime  meridian. 
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Figure  31.  Azimuthal  component  of  the  electrodynamic  winds  in  meters/second 
at  300  km  for  June  solstice  and  0°  magnetic  longitude  as  a  function  of  magnetic 
latitude  and  universal  time.  Contours  intervals  are  20  meters /second.  Eastward 
winds  are  defined  to  be  positive. 
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Figure  32.  Azimuthal  component  of  the  electrodynamic  winds  in 
meters/second  at  1000  km  for  June  solstice  and  0°  magnetic  longitude  as  a 
function  of  magnetic  latitude  and  universal  time.  Contours  intervals  are  20 
meters/second.  Eastward  winds  are  defined  to  be  positive. 
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Appendix  D:  A  Comparison  of  the  IRI, 
Crtiu,  and  SLIM  Models 


A  prime  motivation  for  the  development  of  this  model  for  the  low-latitude,  top¬ 
side  ionosphere  was  to  improve  current  ionospheric  modelling  capabilities  within 
the  U.S.  Air  Force’s  Air  Weather  Service  (AWS).  Since  the  Semi-Empirical  Low- 
Latitude  Ionospheric  Model  (SLIM)  [Anderson  et  al.,  1987]  was  being  consid¬ 
ered  for  incorporation  into  the  Ionospheric  Conductivity  and  Electron  Density 
(ICED)  model  [Kroehl,  1987]  and  ICED  is  the  primary  operational  model  for 
specifying  ionospheric  parameters  within  the  AWS,  the  SLIM  approach  was  re¬ 
viewed  and  its  results  were  compared  to  other  readily  available  empirical  models 
and  known  ionospheric  conditions. 

For  this  comparison,  two  dip  latitudes  (  0°  and  14°N)and  four  local  times 
(00,  06,  13,  and  19  LT)  were  chosen  as  test  cases.  In  the  case  of  the  IRI  and 
Chiu  models,  the  conversion  subroutine  included  in  IRI  was  used  to  convert 
these  dip  latitudes  into  geographic  coordinates.  Since  there  are  no  longitudinal 
variations  in  the  SLIM  [Anderson  et  al.,  1987]  or  Chiu  [Chiu,  1975]  models, 
the  choice  of  255°E  was  arbitrary.  In  addition,  equinox  and  a  sunspot  number 
of  100  (moderate  solar  activity  for  SLIM)  were  used  in  the  model  calculations. 
Figures  33  -  36  show  the  results  of  this  comparison.  Each  figure  shows  the 
electron  density  profile  from  200  km  to  1000  km  as  a  function  of  local  time 
and  dip  latitude.  While  a  number  of  comparisons  could  be  made  using  these 
figures,  the  one  of  primary  interest  was  the  radically  different  trends  between 
IRI  and  SLIM  that  occur  around  800  km.  The  SLIM  profile  shows  an  apparent 
exponential  decrease  in  electron  density  with  height  above  the  peak  in  the  F 
region;  whereas,  IRI  shows  an  apparent  exponential  decrease  and  then  a  much 
less  rapid  decrease  in  electron  density  with  height  starting  between  600  -  800  km. 
IRI  has  been  found  to  underestimat*'  daytime,  total  electron  content  (TEC)  at 
low  latitudes  by  as  much  as  50%,  and  SLIM,  becaiise  it  includes  modelling  of  the 
vertical  E  x  B  drift,  has  produced  significant  improvements  in  TEC  prediction. 
However,  the  SLIM  model  only  computes  O'*'  densities  [Anderson  et  al.,  1987], 
and  H'^  can  become  a  significant,  if  not  the  major,  ion  species  by  1000  km. 
Since  the  AWS  has  requirements  for  modelling  of  the  electron  density  at  least 
up  to  1000  km,  development  of  a  model  which  includes  both  ion  species  was 
undertaken  to  further  improve  ionospheric  prediction  for  altitudes  between  500 
and  1000  km. 
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Figure  33.  A  comparison  of  electron  density  profiles  produced  by  the  SLIM, 
Chiu,  and  IRI  empirical  models  for  0°  dip  latitude  and  324°  magnetic  longitude 
for  equinox  and  moderate  solar  activity,  (a)  0000  LT  (b)  0600  LT  (c)  1300  LT 
(d)  1900  LT. 
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Figure  34.  A  comparison  of  electron  density  profiles  produced  by  the  SLIM, 
Chiu,  and  IRI  empirical  models  for  14°  dip  latitude  and  324°  magnetic  longitude 
for  equinox  and  moderate  solar  activity,  (a)  0000  LT  (b)  0600  LT  (c)  1300  LT 
(d)  1900  LT. 
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Aircraft  data  from  two  flights  during  June  1985,  over  Oklahoma  and 
Arkansas,  are  analyzed  to  explain  the  observed  concentrations  of 
tropospheric  ozone.  Concurrent  vertical  profiles  and  time  series  of 
ozone  with  carbon  monoxide,  odd  nitrogen  and  moisture  are  presented. 
Unique  features  on  this  day  included  pollutant  concentrations  typical 
of  the  planetary  boundary  layer  extending  to  nearly  5  km  altitude,  a 
highly  disturbed  troposphere  over  a  large  area  (about  500  km) 
immediately  following  passage  of  a  squall  line,  and  a  sharp  100  ppb 
gradient  in  ozone  in  the  upper  troposphere  remote  from  convection. 
Correlations  between  ozone  and  other  trace  gases  are  studied  to 
determine  if  the  elevated  ozone  came  from  stratospheric  injection  or 
tropospheric  production.  A  simple  photochemical  model  finds 
significant  tropospheric  ozone  production  possible  under  the 
conditions  observed,  suggesting  that  in  summer  over  the  Central  US, 
tropospheric  ozone  production  greatly  contributes  to  the  observed 
ozone  concentrations. 
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